L  I  B  R.AR.Y 

OF  THE 

UNIVERSITY 

Of    ILLINOIS 

621365 
U655-te 
no.  3 1  ~3G 
cop.  < 


The  person  charging  this  material  is  re- 
sponsible for  its  return  to  the  library  from 
which  it  was  withdrawn  on  or  before  the 
Latest  Date  stamped  below. 

Theft,  mutilation,  and  underlining  of  books 
are  reasons  for  disciplinary  action  and  may 
result  in  dismissal  from  the  University. 

UNIVERSITY    OF     ILLINOIS    LIBRARY    AT    URBANA-CHAMPAIGN 


APR  1  8  11)76 
MAR  2  8 1976 


L161  — O-1096 


Digitized  by  the  Internet  Archive 
in  2013 


http://archive.org/details/dielectriccoated34week 


ANTENNA  LABORATORY 


Technical  Report  No.  J)h 


DIELECTRIC  COATED  SPHEROIDAL  RADIATORS 


Walter  L.  Weeks* 


12  September  1958 


Contract  AF33(6l6)-5220 
Project  No.  6 ( 7-4600 )  Task  *K)572 


Sponsored  "by: 
WRIGHT  AIR  DEVELOPMENT  CENTER 


Electrical  Engineering;  Research  Laboratory 

Engineering  Experiment  Station 

University  of  Illinois 

Urbana,  Illinois 

^Submitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of  Doctor 
of  Philosophy  in  Electrical  Engineering  at  the  University  of  Illinois,  1958° 


>  ACKNOWLEDGMENT 

The  author  wishes  to  thank  all  of  the  members  of  the  Antenna 
Laboratory  staff  for  their  help  and  encouragement,  especially  his 
advisor,  Professor  E.C.  Jordan,  P.E.  Mayes  and  R.L.  Carrel.  Mr.  T. 
Lahey  also  deserves  special  mention  for  his  help  in  the  preparation 
of  the  codes  for  the  ILLIAC.  Financial  support  for  the  investigations 
was  provided  by  the  Wright  Air  Development  Center,  under  Contract 
AF  5220,  for  which  the  author  is  very  grateful. 


ABSTRACT 


Analytical  methods  are  presented  for  the  determination  of  the 
changes  in  the  radiation  patterns  brought  about  by  dielectric  coatings 
on  spheroidal  radiators.  A  description  is  given  of  a  machine  computation 
system  which  makes  the  analytical  results  useful.   This  system  includes 
a  program  for  the  semiautomatic  calculation  of  the  radiation  patterns  of 
spherical  radiators,  a  practical  method  for  the  determination  of  the 
eigenvalues  of  all  of  the  orders  of  the  spheroidal  functions  which  are 
needed  to  obtain  convergent  representations,  and  a  new  method  for  the  cal- 
culation of  the  spheroidal  "radial"  functions.   Numerical  results  are 
included  for  spheres  and  prolate  spheroids,  having  sizes  of  the  order  of 
several  wavelengths,  which  are  excited  by  axially  symmetric  slots. 

The  effects  of  the  dielectric  coating  are  most  marked  when  the 
excitation  is  latitudinally  asymmetric.   One  typical  effect  is  an  increase 
in  radiation  in  the  general  direction  of  the  poles.  For  spherical  rad- 
iators, the  important  effects  show  a  complicated  dependence  on  the  size  of 
the  radiator  and  the  value  of  the  permittivity.   Comparison  of  the  patterns 
of  the  coated  prolate  spheroidal  radiators  herein  calculated  with  the 
patterns  of  spherical  radiators  indicates  that  the  effects  are  similar  in 
both  cases. 
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°  arbitrary  slot  excitation;  used  as  is   0-n^e.     described  above 

(^  Expansion  coefficient  for  spheroidal  angular  functions 

k    ^   Constant  in  dielectric  region  field  representation  in  coated 
°  sphere  problem,  arbitrary  slot  excitation 

|   As  subscript,  means  even  and  odd.  When  both  appear  on  a 
quantity  in  an  equation,  this  indicates  that  the  equation 
actually  represents  two  equations,  one  for  the  even  solutions 
and  one  for  the  odd  solutions 

£e  Vector  mode  function    =  VT«  see  pg.  6k 

e^t.  Vector  mode  function    =  VTe  x  ^      see  #g"  ^ 

\fim%, Constants  in  the  distant  field  expansions  in  the  problem  of 
q,i1ve  the  coated  sphere  with  arbitrary  slot  excitation 


LIST  OF  SYMBOLS  (continued) 

iu^  Scale  factor  corresponding  to  the  coordinate   7A 

Jl_  Scale  factor  corresponding  to  the  coordinate  V 

X4.  Scale  factor  corresponding  to  the  coordinate  d> 

j  Imaginary  unit 

.  Spherical  Bessel  function  of  first  kind  as  used  by  Stratton, 

3U  and  Morse  and  Feshbach 

A  Eigenvalue  symbol 

JL  Eigenvalue  symbol  in  spherical  coordinates 


n. 


Spherical  Bessel  function  of  second  kind  as  used  by  Stratton, 
and  Morse  and  Feshbach 


X.  Semi-focal  distance  in  prolate  spheroidal  coordinates 

■ft  Radius  of  conducting  sphere 

tj_  Outside  radius  of  spherical  shell  of  dielectric 

A 

T  Unit  vector  in  r  direction  in  spherical  coordinates 

tx  Coordinate  symbol.   In  prolate  spheroidal  coordinates,  it  is 
the  reciprocal  of  the  eccentricity  of  the  generating  ellipse. 

IX,  11       -surface  which  coincides  with  the  conducting  spheroid 

11. 


IX       -surface  which  coincides  with  the  outside  of  the 
spheroidal  shell  of  dielectric 


y        Coordinate  surface  symbol „   In  prolate  spheroidal  coordinates, 
it  is  the  reciprocal  of  the  eccentricity  of  the  generating 
hyperboloid 

2         Any  of  the  spherical  Bessel  functions  as  used  by  Stratton,  and 
'L   Morse  and  Feshbach 

cX.  Excitation  coefficient  in  spherical  problems 

cCl  Excitation  coefficient  in  spheroid  problems 

£  Phase  constant,  ^  or  Uj(i<£)3' 

B,        P        for  dielectric  region 

6.         (3       for  free  space  region 


LIST  OF  SYMBOLS  (continued) 

\j,  Kronecker  delta;  equals  one  if    L  =  A  ,    zero  otherwise 

£.  Permittivity  or  dielectric  constant 
B^  E    for  dielectric  region 

£„      £.    for  free  space 

"X  Wavelength 

aj  Permeability 

XJI  Constant  in  the  field  expansion  for  uncoated  prolate  spheroid 

UJ  Frequency  times  att 


Note:  In  section  3,  the  quantities  associated  with  TM  fields  are  denoted 
by  single  primes;  the  quantities  associated  with  TE  fields  are  denoted  by 
double  primes, 


mimunii 


See  note  page  100 


1.   INTRODUCTION 

The  object  of  this  work  is  to  provide  at  least  a  partial  answer 
to  the  following  question;   What  changes  in  the  radiation  pattern 
of  a  moderately  large  antenna  are  brought  about  in  the  event  that 
for  some  reason,  either  within  or  beyond  the  control  of  the  designer, 
the  structure  is  given  a  coating  of  dielectric  material?  The  specific 
antennas  considered  are  conducting  bodies  of  the  spheroid  family, 
excited  by  slots  or  gaps.  In  the  pages  which  follow, an  analytical 
method  is  presented  for  the  determination  of  the  changes  in  the 
radiation  patterns  brought  about  by  spheroidal  dielectric  casings  on 
spheroidal  radiators.  The  basic  geometry  is  indicated  in  Fig.  1. 


Slots  for 
Excitation 


Conducting 
Spheroid 


Dielectric 
Coating 


FIGURE  1.   BASIC  ANTENNA 
It  is  shown  herein  that  numerical  results  can  be  obtained,  in  some  cases 

rather  easily,  with  the  aid  of  a  digital  computer,  and  some  numerical 

and  graphical  results  so  obtained  are  included. 


The  impetus  for  the  study  was  provided  by  the  fact  that  the  shapes 
included  in  the  spheroid  family  provide  reasonable  approximations  to 
certain  space  vehicles  or  missiles,  some  of  which  may  be  given  coatings 
of  (nonconducting)  refractory  material  in  order  to  withstand  the 
temperatures  developed  by  tremendous  speeds  through  the  atmosphere. 
But  this  is  by  no  means  the  only  interesting  aspect  to  the  problem. 
The  effects  of  dielectric  coatings  on  cylinders  and  planes  have  been 
the  subject  of  much  study  in  recent  years  in  connection  with  the 
propagation  and  launching  of  surface  waves,  and  there  is  the  interesting 
question  of  a  possible  connection  or  relationship  between  the  effects 
in  spheroidal  geometries  and  those  in  plane  and  cylindrical  geometries. 
Moreover,  the  general  problem  is  of  interest  to  the  antenna  engineer, 
since  the  effects  of  dielectric  coatings  on  the  one  hand  introduce 
new  possibilities  for  the  control  of  radiation  and  on  the  other  hand 
may  present  difficulties  to  him  when  the  coatings  (such  as  ice  or 
radomes)  are  beyond  his  control. 

The  spheroidal  radiators  are  particularly  attractive  for  a  study 
of  this  kind  since  they  are  the  only  antennas  of  finite  dimensions 
for  which  exact  analytical  formulations  of  the  electromagnetic  boundary 
value  problems  are  possible  (at  least  at  the  present  state  of  knowledge). 
In  spite  of  this,  there  has  been  suprisingly  little  published  work 
concerning  dielectric  coated  spheroids .   The  problem  of  the  scattering 
\   of  a  plane  wave  by  a  dielectric  coated  sphere  having  a  size  near  that 

of  the  first  resonance  has  been  studied  by  Scharfman  and  Aden  and 
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Kerker.  Elliot  has  discussed  the  possibility  of  surface  waves  on 


a  dielectric  coated  hemisphere  set  on  a  perfectly  absorbing  pad. 

As  far  as  is  known  to  the  author,  there  are  no  previous  publications  on 

the  subject  of  the  dielectric  coated  prolate  (or  oblate)  spheroids. 

The  method  employed  in  the  present  work  is  the  classical  method 
for  the  solution  of  electromagnetic  boundary  value  problems .  The 
fundamental  idea  is  that  an  electromagnetic  field  is  uniquely  determined 
in  a  source  free  interior  region  by  the  values  of  the  tangential 
components  of  electric  (or  magnetic)  fields  on  the  bounding  surfaces 
of  that  region*  In  the  antenna  problem,  the  space  exterior  to  the 
antenna  is  effectively  converted  into  an  interior  region  by  the  device 
of  regarding  space  as  bounded  at  one  extreme  by  the  antenna  surface 
and  at  the  other  extreme  by  a  large  enclosing  surface  located  at  a 
distance  from  the  antenna  which  is  so  great  that  the  fields  can  safely 
be  assumed  to  be  those  of  outward  traveling  waves °  In  the  present 
problem,  this  exterior  (rendered  interior)  space  is  not  homogeneous 
(viz.,  it  is  partly  filled  with  dielectric  material).  As  a  result, 
since  solutions  are  easily  obtained  only  in  homogeneous  regions,  the 
problem  is  handled  by  splitting  the  complete  region  into  two  homogeneous 
regions  and  obtaining  solutions  for  each  of  these  regions  by  requiring 
that  the  tangential  components  of  the  electric  and  magnetic,  fields  be 
continuous  at  the  boundary  between  the  homogeneous  regions*   In  outline 
the  solution  proceeds  as  follows:  Maxwell's  equations  for  the  electro- 
magnetic field  are  first  postulated.,  These  vector  partial  differential 
equations  are  transformed  to  ordinary  scalar  differential  equations  by 
introducing  suitable  potential  functions  and/or  by  taking  advantage  of 


known  symmetries.   The  ordinary  differential  equations  so  obtained 
are  of  the  Sturm-Li ouvllle  type,  and  as  such,  with  the  appropriate 
boundary  conditions,  define  characteristic  functions  and  characteristic 
numbers  (eigenvalues)  in  terms  of  which  the  solutions  can  be  expressed. 
The  fields  are  expressed  in  terms  of  these  functions  multiplied  by  (at 
first)  unknown  constants.  Fortunately,  in  the  boundary  value  problems, 
the  number  of  independent  boundary  conditions  which  can  be  applied  is 
exactly  equal  to  the  number  of  linear  equations  which  is  required  for 
the  determination  of  the  unknown  constants.  Unfortunately,  however, 
the  number  of  equations  which  must  be  solved  is  seldom  small. 

In  the  sections  which  follow,  the  foregoing  outline  is  developed 
in  detail.  The  organization  is  such  that  the  problems  are  separated 
first  according  to  whether  or  not  the  fields  have  rotational  symmetry, 
and  then  into  the  specific  coordinate  systems.  In  Chapter  2,  the 
Maxwell  field  equations  in  a  general  coordinate  system  having  rota- 
tional symmetry  are  considered,  and  it  is  shown  that  the  partial 
differential  equations  can  be  separated  in  spheroidal  coordinate  systems. 
The  conditions  for  separability  are  exhibited.  Following  this,  general 
expressions  for  the  field  components  in  such  systems  are  developed  and, 
in  Section  2.1,  these  expressions  are  applied  to  spheres  with  rotation- 
ally  symmetric  excitation  (latitudinal  slots).  Section  2.1.1  includes 
a  review  of  the  solution  for  the  radiated  field  of  a  spherical  antenna 
in  order  to  standardize  the  notation.  The  solutions  for  a  dielectric 
coated  spherical  radiator  are  presented  in  Section  2.1.2.  Also  in 
that  section,  an  expression  is  developed  which  makes  it  possible  to 


evaluate  the  effect  of  a  dielectric  coating,  independent  of  the  type 
of  (rotationally  symmetric)  excitation.   It  is  interesting  to  note  that 
because  of  the  ease  with  which  the  different  orders  of  the  characteristic 
functions  of  spherical  coordinates  (spherical  Bessel  functions  and 
associated  Legendre  functions)  can  be  generated,  the  sphere  calculation 
is  almost  ideally  suited  for  a  digital  computer*  Almost  all  of  the 
results  included  herein  were  obtained  with  the  aid  of  the  University 
of  Illinois  Digital  Computer  (ILLIAC).  The  system  of  calculation  is 
discussed  in  Section  2.1.3. 

In  Section  2.2,  the  general  equations  which  are  developed  in  the 
first  part  of  Chapter  2  are  applied  to  the  prolate  spheroid.  The  problem 
of  the  uncoated  prolate  spheroid  is  reviewed  for  the  purpose  of  standard- 
izing the  notation.  In  Section  2.2.2,  the  solutions  for  the  dielectric 
coated  prolate  spheroid  are  presented.  A  description  of  a  way  in  which 
the  spheroid  calculations  can  be  programmed  for  a  digital  computer  is 
included  in  Section  2.2, 3*  together  with  the  radiation  patterns  for  a 
few  specific  spheroids.  Most  of  the  results  are  new  in  that  they  are 
based  on  function  values  which  heretofore  had  not  been  tabulated .  The 
system  of  calculation  for  the  spheroidal  functions  is  different  from 
that  which  has  been  reported  before  and  it  appears  to  have  some  advantages. 
Chapter  3  is  devoted  to  the  development  of  the  equations  for  spherical 
'  antennas  with  arbitrary  (known)  excitation  —  that  is,  the  restriction  of 
rotational  symmetry  is  removed.  The  equations  for  both  coated  and  un- 
coated spheres  are  developed,  and  it  is  shown  that  the  ratio  of  the 
coefficients  in  the  series  expansions  for  the  TM  fields  is  the  same  as 
that  developed  in  Section  2  for  the  case  of  rotationally  symmetric  excitation, 


The  results  included  in  the  following  pages  show  that  the  effect 
of  a  dielectric  coating  on  a  spherical  radiator  follows  a  rather  typical 
pattern,  the  net  result  usually  "being  an  enhancement  of  the  relative 
magnitudes  of  the  higher  order  Legendre  polynomials  in  the  representations. 
One  effect  then  is  the  introduction  into  the  patterns  of  coated 
spherical  antenna  of  a  given  size  some  of  the  characteristics  of  the 
patterns  of  an  uncoated  antenna  of  somewhat  larger  size,  but  there 
are  other  more  interesting  effects.  The  solutions  for  the  prolate 
spheroid  are  somewhat  less  satisfactory  in  that  the  calculations  require 
a  relatively  large  amount  of  computer  time,  but  more  than  this,  they 
seem  to  be  too  involved  to  provide  insight  into  the  physical  processes 
involved.  The  results  ind: 
that  of  spherical  antennas. 

We  now  turn  to  the  task  of  developing  the  equations  which  apply 
to  general  orthogonal  coordinate  systems  in  the  special  case  that  the 
fields  have  rotational  symmetry. 


2.   SPHEROIDS  WITH  AXIALLY  SYMMETRIC  EXCITATION 


In  the  special  case  that  the  sources  of  the  field  (and  therefore  the 

field  itself)  have  no  variation  with  the  coordinate  (fi       (see  Fig.  2), 

k 

Maxwell's  equations  expressed  in  rotationally  symmetric  coordinate  systems 

(V'>'ir)4>    )  reduce  to  (assuming  e     time  variation)  the  equations 


■a-u. 


—  a  uJA\  Jw  Xv  MA 


(1) 
(2) 

(3) 


together  with  a  similar  independent  set  which  can  he  written  down  "by  inter- 
changing U   and  r I   while  exchanging  the  permittivity  and  permeahility 


FIGURE  2  COORDINATE  SYSTEM 


TABLE  1.   SCALE  FACTORS  FOR  SPHEROIDAL  COORDINATE  SYSTEMS  (u,v,(j>) 


Scale  factors 
and 
related 
quantities 


Coordinate  System 

spherical       prolate  oblate 

(r,(-cos  9),  (j))       spheroidal         spheroidal 


2    2  1/2 

SL    (— ■* )         i  (- 

u  -  1 


2    2  1/2 


u2  +  l 


d-v2)1^ 


2    2 

1  -  V 


2    2  1/2 

1  -  V 


'(1  -  Y2)l/2      I  (  (U2-1)(1-V2)  )X/2  |  (  (U2+1)(1-V2)  )^ 


h*  h 

h 
u 

9     V 

h  h  h. 
u  v  <J> 

fx(u) 

f2(v) 

1-v 


i(l-v2) 

i(l-v2) 

1 

1 

i  (u2-D 

i(u2+l) 

3/    2     2A 
i    (u   -v  ) 

iW) 

3  2 
i    u 

A2 

3  2 

-rv 

3  2 
i   v 
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symbols  according  to  the  rule  &  ——  M         .   The  lk>L   's  are  scale  factors 

of  the  coordinate  system  and  EA.   and  H^   represent  the  components  of 

electric  and  magnetic  fields  in  the  coordinate  system   (."U-.'ir.,  (fr)   . 

If  a  new  variable,    r\=^%v\t    is  introduced,  the  Eqs.  1,  2,  and  5 

can  be  combined  into  the  equation 

1-u.   ^<J>^  "d**-    3V  ^^v  "Bit  X<|> 

Inspection  of  the  table  of  scale  factors  (Table  I)  discloses  that  in 

-n'Tr 

all  of  the  spheroidal  coordinate  systems,  the  quantity  (  ~n — ^ —  )  is 
independent  of  1L    ,  while  the  quantity  (   .  -r  )  is  independent  of  if     , 
so  that  Eq.  h   can  be  written 

^  'S-tA.3-    ^v  Sir* 

Assuming  that  A  has  the  separable  form  A  =  U^OVCir)  ?  it  is 
possible  to  write  the  equation  in  the  form 

V    U        inr   V  (5) 

As  before,  a  glance  at  Table  I  shows  that  the  scale  factors  depend  on  the 
coordinates  in  such  a  way  that   Jk^X^Jkfy  —  -J-, Cu")  ■+-  ^Cv)   (and  as  before 
(  JtA*  )  is  a  function  of  ^l  only  while   (  j^j^  )  is  a  function  of  V 
only)  so  that  the  partial  differential  equation  separates  into  the  two 
ordinary  differential  equations : 
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^U  U     +    (j/^uj^^a  -k)\J    =  o 


V   =  o 


(6) 


(7) 


where  -&   is  the  separation  constant.   Because  of  the  symmetry  of  the 
source  distributions  which  have  been  postulated,  the  v~  —    and  cj)  -  com- 
ponents of  the  fields  must  be  zero  on  the  z-axis.  For  spheroids,  these 
field  components  are  therefore  zero  at  if   —  ±  1    ,   independent  of  tjl 
It  follows  from  the  identity 

with  Eq.  2,  that  these  field  components  involve  the  V  -functions  directly, 
and  therefore  the  following  conditions  on  the  functions  are  implied : 

K  =  0  ,  v  =  ±  1   . 

Equation  7  with  these  boundary  conditions  constitutes  a  Sturm-Liouville 
system,  for  which  it  is  known    that  the  solutions  can  be  represented  by 
characteristic  functions  having  eigenvalues  M.-h    which  give  non-zero 
solutions  of  Eq.  7-   The  general  solution  to  Eq.  6  can  be  written 


Lit ,  and  (Ai 


where   Uo^,  and   (A^  are  the  two  linearly  independent  solutions.   It 
follows  therefore  that  any  solutions  to  the  field  equations  which  meet  the 
specified  conditions  of  symmetry  and  continuity  can  be  written  as  follows: 
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The  electric  field  components  can  "be  found  by  substituting  Eq.  8  in  Eqs. 
1  and  2.   The  result  is 


(9) 


Ev  =  - — l       ILcq^VIU,   -t-^^&oVfe  .    (io- 


Having  thus  derived  the  expressions  for  the  field  components  in  the 
general  spheroidal  system,  we  can  now  proceed  to  the  formulation  of  the 
boundary  value  problems  for  particular  classes  of  spheroids. 

2.1  Symmetrically  Excited  Spheres 

Among  the  spheroids,  the  sphere  is  usually  the  easiest  to  study  be- 
cause of  the  familiarity  of  the  functions  which  are  involved.   This  study 
will  therefore  begin  with  the  sphere  and,  in  fact,  will  present  more  nu- 
merical details  for  it  than  f or vthe  other  spheroids.   The  coordinate 
system  selected  is  that  one  in  which  U.  —  ~t~ ,  v-=—  cos  n   .  i.e.  ( nT,  —  c/>s©>4 

This  selection,  in  place  of  the  usual  (  -f,  G><p  )  system,  makes  the 
sphere  a  more  typical  member  of  the  family  of  spheroids  (the  negative 
sign  in  u  =.—  cos  6  makes  a  right  handed  coordinate  system). 


(r,-  cosG^ct) 
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FIGURE  5  SPHERICAL  COORDINATE  SYSTEM 

With  the  scale  factors  for  this  coordinate  system  taken  from  Tahle  I, 
the  differential  Eqs.  6  and  7  become,  after  replacing   [_J  (ru/)  by   K^  > 


i>R' 


c-^u/Ve  -1S)R    =   0 


(ii) 


and 


i-v*)V"  +^V     =o     . 


(12) 


Although  it  seems  a  pity  to  do  so,  Eq.  12  will  he  transformed  into  the 
more  complicated  hut  more  familiar  associated  Legendre  equation  by  means  of 
the  substitution    Vs  —  Cl-vx)  \  y   as  follows: 


(1 


-^P"  -avp'  +  (4- 


)  P  =o 


(l—v*) 


The  functions  which  satisfy  the  differential  equation,  plus  the  condition 

V3  '~  D  —  u  ,  are  therefore  the  associated  Legendre  function  of  degree  one, 


with  the  eigenvalues  Jfl<,    equal  to  "YiCtu  +  D. 

With  these  eigenvalues,  Eq.  11  is  recognized  to  be  a  Bessel-type  differ- 

7 
ential  equation  which  has  solutions  of  the  form 

R    =  o,^^^    -+■     A^N^pf) 

(  he   -  tU^+I)  identified  as   Mjg;  -1    ,  so  -f  -  *  ^+i    ,  where  p  is 

a.  "+ 

the  function  order  )  (3  =  U)  >m,£  .   These  functions  can  now  he  inserted  into 

the  general  Eqs.  8  and  10  for  the  field  components  which  are  tangent  to  the 

spherical  surfaces   In  doing  so,  it  is  convenient  to  use  the  notation  of 

Schelkunoff  and  introduce  a  radial  function 

*—v^  X  /L—n.+.L 

3-> 

where    ZL^^^   is  any  one  of  the  three  kinds  of  Bessel  functions.   In  terms 
of  these  functions,  the  equations  for  the  tangential  field  components  become 


and  from  Eq.  2 


where  A 

r- az    . 


Ik 

2.1,1  Boundary  Value  Problem  of  the  Conducting  Sphere  (no  coating) 

Unique  solutions  of  boundary  value  problems  in  regions  exterior  to 
a  spherical  surface  are  obtained  from  the  values  of  tangential  L   (or 

(— |   )  over  the  spherical  surface,  together  with  the  radiation  condition 
at  a  larger  spherical  surface  whose  radius  tends  to  infinity.   Since  the 
radiation  condition  implies  that  the  fields  are  propagating  in  the  di- 
rection of  increasing  -f   ,  we  note  from  Eq.  13  that  the  coefficients 

O-n.  and  Jcrr^   must  be  related  so  that  Jr  l  =  -a  Oj^  ,    in  order  to  give  the 
proper  function  to  represent  the  outward  traveling  waves , 
A  £_  A 

H„<p-rt  =  Jj?^  -  \  Htf*-5  • 

Thus  the  expressions  for  the  fields  exterior  to  a  conducting  sphere  can 
be  written  in  the  form 

^  -r         -yu 

A 

^  =  ~^b  E^H>P>  .  (16) 

The  coefficients  in  these  expansion  are  determined  by  requiring  that  the 
series  expansion  Eq.  l6  represent  the  electric  field  at  the  surface  of  the 
conducting  sphere  (   ■+■  —  -fc  )  *  This  latter  electric  field  is 
of  course  zero    except  over  the  axially  symmetric  slots.   The  field 
in  the  slot  is  presumed  "to  be  known,  and  if  it  is  represented  in  the  form 
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)  ,*_  P«     .  (17) 


the  orthogonality  properties  of  the  Legendre  functions  make  it  possible  to 
obtain  the  following  equation  for  the  coefficients : 


=  jL  l  R.1  a 


SIOt3       ' CJL         f^ 


where 


TcPjAir 


(18) 


X   <F^)a'Air    —      a/yTb+y 


(19) 
2.^  +  1 

The  expansion  Eq.  16  evaluated  at  -f*  •=.  -f;  must  he  identical  to  the  ex- 
pansion,Eq.  17;  this  implies  that  the  coefficients  must  he  equal  term  by 
term.  Thus  the  unknown  coefficients  in  Eq.  16  can  be  found  in  terms  of  the 
excitation  coefficients  oC^,     from  the  equation 

a~=--^^  •  (CO) 

The  formula  for  the  computation  of  the  fields  reads  therefore 

H,  ^-d^-r,  ]T      A1     (f  E^Efrv>Ay  >  H^^V)  (2D 

In  the  calculation  of  the  distant  field,  the  asymptotic  form  for   H.^f3 
can  be  used: 


H  w 


x— *-©° 


e 
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Thus,  the  expression  for  the  calculation  of  the  radiated  field  is 

To  display  the  variation  of  the  field  with  angle   9    (the  pattern), 
only  the  quantities  to  the  right  of  the  summation  sign  are  of  interest 
This  much  of  equation  (22)  was  programmed  for  calculation  by  the  ILLIAC. 
Some  of  the  details  of  the  solution  by  the  digital  computer  are  presented 
in  Section  2.1.3- 

In  a  desk  calculation,  the  "book  Tables  of  Spherical  Bessel  Functions 
can  be  used  to  advantage.  Therein,  the  function  Z^Cx)  =  C-^-T"  zi^^is 
tabulated.  With  these  tables,  the  following  relationships  are  helpful: 


Z^P-t)  =  (5+-  zjtprt 

A 

<£+■ 


2.1.2  Boundary  Value  Problem  of  the  Coated  Conducting  Sphere 

The  equation  for  the  distant  field  intensity  is  basically  the  same 
whether  or  not  the  radiator  is  dielectric  coated.   However,  the  coefficients 
in  the  field  expansions  are  of  course  different  and  these  are  obtained  from 
the  solution  of  the  boundary  value  problem.  One  measure  of  the  effect  of  a 
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FIGURE  5  DIELECTRIC  COATED  SPHERICAL  RADIATOR 
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dielectric  coating  on  the  surface  of  a  sphere  can  be  obtained  by  comparing 
the  coefficients,  as  given  by  Eq.  20  with  the  corresponding  coefficients 
in  the  coated  problem.   This  is  done  in  the  present  section. 

When  the  sphere  is  dielectric  coated,  we  must  recognize  two  regions: 
a  region  (i)  consisting  of  the  spherical  shell  of  dielectric  surrounding 
the  conducting  sphere,  and  a  region  (il)  consisting  of  the  remainder  of 
the  exterior  space.   In  region  (i),  the  Eqs.  13  and  1^  apply  directly.   In 
region  (il),  the  equations  have  the  same  form  as  Eqs.  15  and  l6,  which  can 
be  rewritten  as  follows: 

l-L    =AVcrH'.rvoP>)  (23) 

Again,  the  series  expansion  for   t_       ,  evaluated  at  -T"  =  "tj   ,  must 

be  equal  to  the  series  expansion  for  the  applied  field  at  that  surface,  which 

latter  is  taken  to  be  that  of  Eq.  17.  Equating  the  coefficients  in  these 

two  expansions,  term  by  term,  gives  the  following  relationship  among  the 

coefficients : 

a  a  f 

a^J.,'/.fU>  4-  J^N_,./fv^  =-^£i+,«<tu    .      (25) 

Equating  the  tangential  component  of  the  magnetic  fields  at  the  surface 
■+-  =  ■*-__  ,  by  equating  (13)  evaluated  at   ^   to  (33)  evaluated  at   "f"x   , 
gives  the  following  relationship: 
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a-,. 


A  A  A 


(26) 


The  corresponding  condition  on  the  tangential  component  of  the  electric 
fields,  applied  by  equating  (ik)   and  (2k)   evaluated  at  -t  ■■=.  -|%    , 
gives  the  following  relationship: 

*  A  A 

^  Jj^  +  AcNL^^  =  f  c-  H^P-«    •        (27) 

Equations  25,  26,  and  27  constitute  sets  of  three  equations  in  the  three  sets 
of  unknown  expansion  coefficients.  The  augmented  matrix  for  this  system  of 
equations  is 


A 
A 


-JuJE.-CoC 


N>*o       o 


Thus,  we  find  an  equation  for  the  expansion  coefficients  in  region  (il)  as 
follows ; 


A  A  /\  /\ 
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where 


A^= 


o 

A 


A  A 

A  A  A 


The  quantity  in  traces  in  Eq.  28  is  the  Wronskian  determinant  of 

/\      A 

(   J   IN^ )  •  1^is  quantity  can  he  determined  from  the  asymptotic 
solutions  of  Eq.  11.  The  result  is 

With  this  result,  the  equation  for  the  distant  field  coefficients 
takes  a  simpler  form  as  follows  t 


(29) 


^  =  -4^^^ 


A. 


(30) 


Now  a  measure  of  the  effect  on  the  radiation  pattern  of  the  dielectric 
layer  on  a  conducting  sphere  may  be  found  from  studying  the  ratio  of  the 
coefficients,  (X^_  (Eq.  20),  which  apply  to  the  uncoated  sphere  to  the 
coefficients,  C^,  (Eq.  30),  in  the  situation  that  the  conducting  sphere 
has  the  same  size  and  the  same  excitation  hut  has  a  dielectric  layer.  The 
ratio  is  as  follows ; 


cu.    '""'     i    ftl-Ccp.*,) 


(3D 


^    '  '  **  A^ 


(32) 
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It  is  interesting  to  study  these  ratios  in  limiting  cases.   For 
example,  if  "both  spherical  surfaces  are  large  enough  so  that  the  Bessel 
functions  can  be  replaced  by  their  asymptotic  forms  for  large  arguments,  i.e., 

A 

J^ >  cos  Cp-r-{c*t+xN 

A 

e 


re— »  r 


the  result,  after  a  fair  amount  of  rearrangement,  can  be  put  into  the  follow- 
ing form; 

c^      ~a^    e . ^,-**\ 

a.       °  sivu  p.c-r,.-^,)  —  ^ ^  cos  ftto-r.) 

It  is  noteworthy  that  this  ratio  is  independent  of  n.  This  indicates  that 
for  such  large  spheres,  the  dielectric  coating  will  not  influence  the 
pattern  variation  with  angle,  but  the  approximations  may  not  be  valid  for  the 
highest  orders  needed  in  this  representation. 

On  the  other  hand,  if  (  (3>N~  )  is  very  small,  we  can  take  as  approxi- 
mations the  first  terms  in  the  power  series  representations : 


-r\. 


C  "XT**  -*"  1^  I 


n™  =  - 


C3.Tb)j  V 


-■n. 
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Then,  after  a  fair  amount  of  algebra,  the  ratio  of  coefficients  becomes 


For  such  small  spheres,  only  the  term  having   tl   equal  to  unity  is 
important : 


a,   -L     a. 

C,      3£r 


^F  +  i^1"^1] 


It  should  be  mentioned  here  that  there  is  a  danger  of  making  an  invalid 
observation  based  on  the  latter  equation.  This  equation  applies  in  the 
event  that  the  electric  field  is  the  same  in  both  the  coated  and  uncoated 
state.  Normally  in  practice  the  constant  would  be  the  voltage  across  the 
gap,  in  which  case  the  electric  field  would  be  less  with  the  dielectric 
present. 

For  the  intermediate  (and  more  interesting)  range  of  values  of  (S-r- 
the  ratio  must  be  studied  numerically.  This  was  done  and  the  results  are 
discussed  in  Section  2.1.3. 

2.1.3  Numerical  Computations  Associated  with  the  Spherical  Radiator 

Because  of  the  relative  ease  with  which  different  orders  of  the 
spherical  functions  can  be  generated,  the  sphere  calculation  is  almost 
ideally  suited  to  a  digital  computer.  The  functions  which  are  involved 
(spherical  Bessel  functions  and  associated  Legendre  polynomials)  satisfy 
very  convenient  recurrence  relationships  as  follows : 
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A 


and 


?   CjO  =   3.^+1  7  (x)  "  7  <*-> 

z,>>  ■=  ^  z>>  - z^x  =^^z^-^zr^l 

P1  (%)  -    ±-  \  ca^H-^^P1^)  -(ti  +  i^PV)! 

Thus,  given  the  two  lowest  orders  of  the  functions  (which  involve  only  sines 

and  cosines),  the  higher  orders  are  easily  generated  in  a  digital  computer 

"loop".  The  calculations  described  herein  were  "based  on  these  recurrence 

relationships  -  Qeneration  of  the  functions  was  done  entirely  by  the 

machine. 

With  Legendre  polynomials  and  the  spherical  Bessel  functions  of  the 
A 

second  kind  (  j\J  )  this  upward  recurrence  process  can  apparently  he  carried 
out  to  almost  indefinitely  high  orders  of  the  functions  without  serious  loss 
of  accuracy.  However,  in  the  calculation  of  the  spherical  Bessel  functions 
of  the  first  kind   u^  ,  there  is  a  serious  loss  of  accuracy  if  carried 
out  in  a  straightforward  upward  recurrence  process.   In  fact,  the  figures 
begin  to  lose  significance  for  those  orders  (  ru  )  which  exceed  approximately 
twice  the  argument  --  and  this  is  approximately  the  number  of  terms  re- 
quired for  the  series  expansions  for  the  fields  to  converge  satisfactorily. 
Hence,  an  alternative  system  is  employed  for  these  functions.   The  relation- 
ship, 


T   _  aTl-h  i .  7  _  7 


2k 


is  transformed  into  a  continued  fraction  form  for  use  in  a  downward  scheme 
as  follows : 


^j-A-t-i  —  a.n-t-1  L-li.       _  ^ 


i* 


X, 


J, 

— 7v 

T 


.Jtl     .    ^+1 


I 


1 


a.-n.  +  i 


1 

QjfL+3 

1 

^ 

\ 

\ 

\ 


\ 


3.(TL+m)-l   JC 


J, 


n+m-i 


/A 

Now  if  -w\.  is  taken  sufficiently  large,  the  ratio  J^+^w/  J^+Yrv. 
negligihle  in  comparison  to  the  term  a,  (n+m -!)/;&   ("because  of  the  nature 
of  the  functions  of  high  order  and  relatively  small  argument),  and  hence  a 
good  approximation  for  this  ratio  is  zero.  With  this  approximation,  the 


ratio 


I 


is  first  calculated  and  then  employed  to  cal- 


culate  the  lower  order  ratio.  This  process  is  continued  so  as  to  calculate 
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and  store  the  ratios  of  successive  orders  of  functions  from  some  high  order 
to  some  low  order  (say  n  -  1    )  for  which  the  value  of  the  function  in  the 

A 

denominator  is  readily  obtained  (  JT(v.)  —  sin.%   ).  The  functions  them- 
selves are  obtained  from  the  ratios  by  successive  multiplications.  When 
generated  in  this  way,  all  except  the  highest  several  orders  calculated  can 
have  full  machine  accuracy.  However,  the  method  breaks  down  if  the  argu- 
ment coincides  exactly  with  a  zero  in  one  of  the  (lower)  orders.  Thus,  it 
is  safest  to  use  the  upward  recurrence  form  in  the  oscillating  region  of  the 
functions  and  the  continued  fractions  scheme  in  the  monotone  region.  An 
alternative  is  the  use  of  the  recurrence  form  in  the  downward  sense.  This 
latter  method  is  satisfactory  except  for  the  arguments  which  correspond  to 
the  zeros  of  the  lowest  order  function  computed  (integral  mutliples  of  2jt 

A 

if  the  process  is  carried  to  Joto  )• 

Except  for  this  difficulty  with  the   J-a,  functions,  a  completely  auto- 
matic digital  computer  program  could  easily  be  written,  in  that  the  decision 
as  to  the  number  of  terms  to  be  computed  could  be  made  by  the  machine  (by  a 
convergence  test).  The  actual  program,  described  below,  which  was  written 
for  the  ILLIAC  is  semi-automatic  in  the  sense  that  the  programmer  makes  a 
guess  concerning  the  number  of  terms  required  for  convergence,  and  on  this 

A 

basis,  a  sufficient  number  of  J^  functions  is  made  available.   If  the 
guess  which  is  made  is  too  low,  the  program  stops  on  a  special  stop  order. 
(The  number  of  terms  required  appears  to  be  of  the  order  of  twice  the  argu- 
ment of  the  Bessel  functions  (  fi"f  )  or  perhaps  ten  times  the  radius  in 
wavelengths.  Actually,  orders  of  the  functions  as  high  as  sixteen  times  the 
radius  in  wavelengths  were  computed  in  the  program,  and  this  number  never 
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failed  to  provide  a  safe  margin). 

In  the  calculation  of  the  radiation  patterns,  the  constants  and  the 
■^  variation  in  Eqs.  22  and  23  are  suppressed  and  an  expression  having 
the  form 


y-  =  yl  r  kp: 


is  evaluated.  The  quantity  A^  is  of  course  a  complex  number  which  is 
evaluated  as  follows,  in  the  uncoated  and  coated  problems  respectively: 

A  A 

JT  +  re  " 


where 


and 


in  which 
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FIGURE  6  SIMPLIFIED  FLOW  DIAGRAM  FOR  RADIATION  PATTERN  CALCULATIONS 
(SPHERICAL  RADIATOR) 
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A  simplified  flow  diagram  for  the  machine  calculation  of  the  radiation 
patterns  is  presented  in  Fig.  6.  The  computer  program  is  written  in  such  a 
way  that  for  a  given  angle  in  the  field,  additional  terms  in  the  summation 
(on  tl  )  are  obtained  in  groups  of  four.  This  summation  process  is  carried 
on  to  higher  and  higher  orders  until  a  comparison  with  the  cumulative  sum 
shows  that  the  last  four  terms  computed  is  satisfactorily  small  in  com- 
parison to  the  cumulative  sum.  The  excitation  coefficients  (  oC^  )  are 
obtained  by  a  five  point  Simpson's  rule  integration  process  on  the  assumption 
that  the  electric  field  is  constant  over  a  slot  whose  width  subtends  two 
degrees  of  latitude.   The  machine  time  required  (using  a  floating  point 
interpretive  routine)  for  a  single  pattern  calculation  involving  thirty- 
five  points  in  the  field  is  about  five  minutes  (for  greater  accuracy  than 
is  required  for  a  cathode  ray  tube  display).  If  several  patterns  are 
computed  by  varying  the  parameters  while  the  program  is  still  in  the  machine, 
the  average  time  per  pattern  is  less. 

Some  of  the  radiation  patterns  which  were  so  calculated  are  displayed 
in  Figs.  13  to  20.  The  radiation  patterns  of  a  given  spherical  radiator 
with  and  without  a  dielectric  coating  are  presented  side  by  side  for  con- 
trast. 

The  number  of  radiation  patterns  included  herein  is  not  large,  since 
it  appears  to  be  ultimately  more  useful  and  more  economical  to  present 
certain  key  data  pertaining  to  spherical  radiators,  from  which  additional 
c&lculations  can  be  made  with  a  minimum  of  effort.   Thus,  Table  k   gives  the 
excitation  coefficients  (calculated  as  mentioned  above  in  the  hope  that 
they  will  be  useful  in  future  impedance  calculations)  for  slots  whose  posi- 
tion varies  from  ten  to  ninety  degrees  latitude.  Another  quantity  of 
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interest  is   |  rl^fvnl  }   for  the  uncoated  radiator.   The  square  of  this 

quantity  is  plotted  as  a  function  of  ~^/\     for  a  given  Tl   in  Fig.  28:. 

The  maxima  in  these  curves  show  the  sizes  for  which  the  sphere  becomes 

qua si -re sonant  for  the  various  orders  of  the  functions.   The  same  quantity 

as  a  function  of  "Tx.  ,  for  different  values  of  P0^   is  given  in  Figs. 29 

A   /  I  A  i    i  ^.     A '   I  ^Vl a 
to  39,  and  in  Table  5  with  the  quantities   <X\./   I  Hx  |  and  N_/ 1 nTJvhich 

are  needed  in  any  pattern  calculations.   Thus,  with  these  data,  and  a  good 

table  of  Legendre  polynomials,  a  radiation  pattern  calculation  can  be  made 

relatively  easily  with  a  desk  calculator. 

Similarly,  some  key  data  are  included  for  the  dielectric  coated 

spherical  radiator.  As  was  mentioned  earlier,  the  ratio  of  coefficients 

(coated  to  uncoated)  given  by  Eq.  32  is  studied  in  some  detail.  However, 

it  appears  to  be  ultimately  more  useful  to  present  the  data  for  the  real 

parts  and  the  imaginary  parts  separately,  since  with  this  done,  given  the 

coefficients  in  the  expansion  for  the  uncoated  case,  one  can  immediately 

find  those  for  the  coated  sphere  from  the  Tables,  Thus,  in  Figs,  1*0  to  63, 

and  Table  6,  the  results  of  a  study  of  the  equations 

Rec,  =  e$  1  fc<w r  JL  . 

are  presented.  These  ratios  were  evaluated  as  a  function  of  the  size  of  the 
conducting  sphere,  the  dielectric  constant,  and  the  thickness  of  the  coating. 


50 
Three  sets  of  calculations  were  carried  out,  in  each  of  which  one  of  the 
above  quantities  was  varied  while  the  other  two  were  held  constant.  Note 
that  there  is  a  rather  characteristic  pattern,  which  is  test  exemplified  by- 
Fig.  b2  +     The  characteristic  peak  which  appears  always  indicates  an  enhance- 
ment of  the  magnitudes  of  the  higher  order  Legendre  polynomials  in  the  field 
expansion.   The  peaks  (orders  of  greatest  enhancement)  grow  higher  and 
move  to  higher  orders  as  either  the  value  of  the  permittivity  or  the  thick- 
ness of  the  coating  is  increased  (for  a  given  size  of  conducting  sphere). 
There  are  other  interesting  features  in  the  graphs,  some  of  which  can  be 
misleading.   It  hardly  needs  to  be  pointed  out  that  the  curves  will  exhibit 

A 

peaks  in  the  vicinity  of  the  zeros  of   /  -,yj   j  ^^   ^e  Pea^s  may  not  be 
significant  as  far  as  the  patterns  are  concerned  since  the  contribution 
of  those  orders  to  the  uncoated  radiation  pattern  may  be  negligible.  The 
peaks  are  generally  higher  in  the  curves  of  Ke  Cv /Ke(X,(X,  but  again  the 
effect  in  the  patterns  may  be  less  marked  if   Re  (X,^     does  not  itself 
make  an  appreciable  contribution  to  the  uncoated  pattern.   It  can  also  be 
seen  that  the  extent  of  the  effect  depends  upon  the  nature  of  the  excitation. 
For  example,  suppose  that  there  is  a  rather  sharp  peak  at  a  particular  odd 
(even)  order  of  the  functions  $   then,  there  will  be  a  marked  effect  if  the 
odd  (even)  orders  are  excited  but  practically  none  if  the  odd  (even)  orders 
are  not  excited.   (Even  orders  are  not  excited  by  slot  fields  which  are 
symmetric  about  the  equator j  odd  orders  are  not  excited  by  excitations  which 
are  anti-symmetric  about  the  equator,  such  as  a  pair  of  slots  having  the 
same  latitude,  one  in  the  northern  hemisphere  and  one  in  the  southern  hemi- 
sphere, excited  one  hundred  and  eighty  degrees  out  of  phase.)   In  general, 
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the  effects  are  greater  when  the  excitation  is  such  that  all  orders  are 
present,  and  one  typical  result  is  an  increase  in  the  radiated  field  in 
the  general  direction  of  the  poles  (z-axis).   This  effect  arises  because 
of  the  enhancement  of  the  higher  orders  of  the  Legendre  polynomials.   If, 
howeve^  certain  of  the  adjacent  higher  order  ratios  have  comparable  mag- 
nitudes but  opposite  signs,  some  surprising  results  may  arise  as  may  be 
seen  in  some  of  the  radiation  patterns. 

We  next  turn  our  attention  from  spherical  radiators  to  prolate 
spheroidal  radiators,  excited  by  axially  symmetric  slots. 

2.2  Symmetrically  Excited  Prolate  Spheroids 

The  -dielectric  clad  spheroid  which  JLs  probably  of  most  practical 
interest  at  the  present  time  is  the  prolate  spheroid.   In  the  prolate 
spheroidal  coordinate  system,  the  values  of  the  coordinate   1A.   are  the 
reciprocals  of  the  eccentricities  of  the  ellipses  which  generate  the  pro- 
late spheroids  (so  1  £  xc  ^  cO  ).  Th.e  orthogonal  -\r  surfaces  are 
hyperboloids  of  two  sheets  and  the  values  of  the  coordinate   or   correspond 
to  the  reciprocals  of  the  eccentricities  of  the  generating  hyperbolas 
(-1  ^  V-  —  1  )   The  relationship  of  this  coordinate  system  to  the  spheri- 
cal system  can  be  illustrated  by  writing   u  —  —  cos  ©  ,   where   0   is  the 
(polar)  angle  which  is  defined  by  the  asymptotic  cones:  thus,  as  1(  — ^*=c> 
so  that  the  surfaces  u  =  constant  become  spherical,  the  angle  0  becomes 
the  usual  polar  angle  of  spherical  coordinates   The  quantity  JL      is  the 
semi-focal  distance. 
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FIGURE  7  PROLATE  SPHEROIDAL  COORDINATE  SYSTEM 
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With  this  identification  of  coordinate  variables,  the  results  of  the 
analysis  presented  at  the  beginning  of  Section  2  are  immediately  appli- 
cable.  It  is  clear  from  the  results  in  Table  I,  and  Eqs .  6  and  7,  that 
the  differential  equations  which  apply  to  the  prolate  spheroidal  system 
are  as  follows : 


Ctjl* 


-dU"  ■+  <<p*uf -J^LI  =■  °        (36) 


The  functions   [_J        and  \/   are  only  incompletely  tabulated.   Hence, 
the  first  task  in  the  spheroid  problem  is  to  devise  a  system  for  the  cal- 
culation of  these  functions,  including  the  determination  of  the  family  of 
eigenvalues,  Jkj^      ,   so  that  the  functions  \/^   comprise  a  complete  ortho- 
gonal set* 

The  system  of  calculation  finally  adopted  was  selected  with  an  eye 
toward  the  most  advantageous  use  of  the  ILLIAC  and  its  library  of  service 
routines.  An  outline  of  the  system  is  presented  in  Section  2.3-3   The 
details,  together  with  some  numerical  values,  will  appear  in  a  separate 
technical  report  . 

Essentially,  the  \J   -functions  are  calculated  from  a  series  of 
associated  Legendre  polynomials : 


Vt  =  ci-^  YL  <«*.*£■  .       (38) 

TV 


where  the     ^^.n  satisfy  the  recursion  relationship: 


-Cfilf 


i 


i^x-a. 


1  !  Ca.w.-*-3Xajr\.-t-5")       ^^  +  SL 


(We  note  here  that  the  eigenvalues  Jk^  depend  on  (3#-  so  that  the  matching 
of  the  field  components  at  a  boundary  between  dielectrics  is  complicated.) 

The  "radial"  functions,   \ J   ,  are  calculated  from  a  series  of  spherical 

Bessel  functions,  the  coefficients  in  which  satisfy  a  recursion  relation 
very  similar  to  that  given  above. 

2.2.1  Conducting  Prolate  Spheroid  without  Dielectric  Coating 

Referring  back  to  Eqs.  8  and  10,  after  inserting  the  scale  factors  for 
the  prolate  spheroidal  coordinate  system,  it  can  be  seen  that  the  field  ex- 
pressions are  as  follows  (  ^IXj^  is  the  solution  of  Eq.  36  which  satisfies  the 
radiation  condition): 

U  =   ' " L   /  ^'UfcsVfc  (to) 

As  was  done  in  the  case  of  the  spherical  radiator,  the  tangential  component 
of  the  electric  field,  F.^  ,  is  specified  to  be  zero  over  the  (conducting) 
'surface  7^  =  ix,  ,  except  over  the  azimuthal  slots.  The  field  on  that  sur- 
face is  then  expanded  in  the  orthogonal  set  of    V   -functions  as  follows: 
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It  follows  from  the  boundary  conditions  and  Eq«  37  that  the  statement  of 
orthogonality  among  the  V  -functions  is  as  follows : 


J_x      (l-vl) 


O         if     A-    *:   TL 


so  that 

E,  c-^-v1^  V.   A- 


0  Js/o 


^   —  _     1 ::' LA-y±± 


^-1        ci-v*) 


(42) 


The  actual  evaluation  of  the  coefficients  is  done  in  terms  of  the  series 
(38)  and  (kO)   so  that 
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In  the  numerical  calculation  of  the  \J   -functions,  there  is  an 
arbitrary  constant  to  be  chosen.   It  appears  from  the  latter  equation 
that  a  very  convenient  choice  would  be  one  which  leads  to  the  following 
equation: 

y  k,  f <  p;^  ^  =  y  <&  ™™*  =  1 

The  arbitrary  constant  is  so  chosen  in  the  calculation  described  in 
Section  2.2.3. 

Also,  in  the  calculation  of  the  distant  field,  it  is  convenient  to 
have  a  familiar  asymptotic  form  for  the  U  -functions  as  (31  u  — *  (3l~-><=> 
The  arbitrary  constant  in  the  Bessel  function  series  for  these  functions 
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FIGURE  8  SPHEROIDAL  RADIATOR 


Dielectric 
Region  (I)  (€1) 


Region  (II)  (€Q) 


FIGURE  9  SPHEROIDAL  RADIATOR  WITH  DIELECTRIC  COATING 
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is  chosen  so  that 

A 

that  is,  the  (A  -function  is  made  asymptotically  equal  to  the  corre- 
sponding order  of  spherical  Bessel  function. 

With  these  normalizations,  the  equation  for  calculating  the  dis- 
tant field  becomes 

2.2.2  Conducting  Prolate  Spheroid  vith  Dielectric  Coating 

Let  the  conducting  spheroid  coincide  with  the  surface  "U  —  Ua  , 
and  let  the  dielectric  coating  extend  from  the  surface  Ik   —  ut  to  the 
surface  "U.   ^^   .  The  dielectric  region  will  he  designated  (i)  and 
the  remainder  of  exterior  space  will  he  designated  (il).  Differential 
equations  of  the  type  (36)  and  (37)  apply  to  both  regions;  however,  the 
propagation  constants  and  the  eigenvalues  are  different  for  the  two 
regions.  This  complicates  the  boundary  value  problem  since  the  eigen- 
f unctions  on  one  side  of  the  boundary  are  different  from  those  on  the 
other  side  of  the  boundary,  which  means  that  the  coefficients  cannot  be 
equated  term  by  term. 

In  region  (i),  the  equations  for  the  tangential  components  are  the 
same  as  (8)  and  (10),  with   £  =  £j.    , 


(*3) 
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E„  =-[c^-^)(i-V)]i7[<  ^  ^  +  J^  AT^^Vt     .      (lop) 

In  region  (II),  the  equations  are  of  the  type  of  (\o)   and  ( I4-1) ;  however, 
in  writing  the  expressions,  the  eigenvalue  symbol,  Jk     ,   will  be  re- 
placed by  the  letter  -j   ,  to  underscore  the  fact  that  the  functions  are 
different  in  the  two  regions : 


o 

^  u;£or        ' '  ^  1^'"   '  (laa) 

The  field  is  presumed  to  be  known  over  the  inner  (conducting)  Spheroidal 
surface,  xl±   ,  (the  tangential  electric  field  is  zero  except  over  the  slots) 
and  is  expanded  in  a  series  of  V  -functions : 


E       =  f^  =  -(c^-^)u-^)]""  V  ^\4 


The  coefficients,  otj^   ,  are  determined  in  the  usual  way  by  a  multiplication 
by   V-  (l-^r1    ,  with  the  result: 


oCA-     d 
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The  actual  computation  of  these  coefficients  can  be  done  with  Legendre 
polynomials  and  (X    -coefficients  as  indicated  in  Section  2.2.1,  under 
Eq.  k2.     Thus,  the  application  of  the  "boundary  condition  at  the  surface 
U!  leads  to  a  set  of  -&   equations  in  the  unknown  coefficients  cl^ 
and   -&I?.   ,  as  follows i 

TZ  Av.  (w° 

Two  more  sets  of  equations  are  obtained  from  the  boundary  conditions 
at  the  surface  lx   =^^  •  The  condition  on  the  magnetic  field  gives  the 
result 


? 


while  the  condition  on  the  electric  field  gives  the  equation 

Equation  1+4  can  be  employed  to  remove  the  variable,  Xfj^     ,  from  Eqs.  k-5 
and  k6  as  follows : 
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The  quantities    V^   and    \A'   can  be  replaced  as  follows: 


The  reason  for  the  latter  step  is  that  now  the  corresponding  orders 
of  the  Legendre  polynomials  can  he  equated  term  by  term,  whereas  this 
is  not  possible  for  the  V  -functions.  When  this  is  done,  and  the 
substitutions 

are  introduced  to  save  writing,  the  result  is 


in 

These  latter  equations,  (kj)   and  (kQ) ,  must  be  solved  for  the  unknown 
coefficients   C^   (and   CX.^  ).  The  equations  are  valid  for  each  value 
of  yv.  ,  so  that  each  of  the  equations  can  be  written  down  the  number  of 
times  which  is  required  to  solve  for  the  Jk  -4-  A      unknown  coefficients.  Of 
course,  the  solutions  to  the  equations  are  strictly  correct  only  for  an 
infinite  set.  However,  the  series  expansion  of  the  fields,  and  those  for 
the  V  -functions,  are  presumed  to  be  convergent  representations.  Hence, 
it  will  be  permissible  in  an  actual  computation  to  replace  the  infinite 
sums  by  sums  having  a  finite  number  of  terms.  Then,  the  number  of  equa- 
tions which  must  be  handled  is,  of  course,  at  least  twice  the  number  of 
terms,  N,  which  are  required  for  a  sufficiently  accurate  distant  field 
representation.  The  fact  that  the  ^^  -coefficients  peak  rather  sharply 
about  the  value,  o-jkJk,      ,  (i.e.  the  \J    -function  of  a  given  order 
closely  resembles  the  Legendre  polynomial  of  the  corresponding  order)  means 
that  a  finite  set  of  equations  can  be  essentially  decoupled  from  the  in- 
finite set,  and  this  finite  set  may  not  be  appreciably  greater  than  the 
minimum  number  required  (2N).  Hence,  in  practice,  one  can  start  with  some 
small  number  of  equations  (say  a  number  of  equal  to  three  times    (31   ), 
obtain  the  solutions,  then  increase  the  number  successively  in  steps  of  four 
and  repeat  until  the  following  two  conditions  are  satisfied:  l)  The  co-, 
efficients  determined  do  not  change  appreicably  as  the  number  of  equations 
is  increased  and  2)  the  magnitude  of  the  highest  order  coefficients  deter- 
mined are  some  small  percentage  of  the  largest  (the  latter  condition  will 
guarantee  the  former,  but  the  largest  coefficient  may  sometimes  be  hard  to 
find). 


k2 

It  is  worthy  of  special  mention  that  the  odd  orders  and  the  even 
orders  are  completely  decoupled,  and  this  good  fortune  cuts  the  number 
of  simultaneous  equations  in  a  set  in  half.   However,  the   O-a.  }s        and 

C,  's  (and  perhaps  the  °Cfc '«  )  are  complex  numbers,  and  the  necessity 

for  equating  the  real  parts  and  the  imaginary  parts  separately  doubles 
the  number  again. 

Once  the  C^  's  have  been  determined,  Eq.  4la  can  be  used  to  deter- 
mine the  distant  field,  wherein  the  U,3 -function  is  replaced  by  its 
asymptotic  form,  as  was  done  in  Eq.  k^. 


2.2.3  Numerical  Calculations  Associated  with  Prolate  Spheroidal 
Radiators 


There  are  two  very  serious  obstacles  to  the  carrying  out  of  the 
numerical  calculations  for  spheroidal  radiators   The  first  is  the 
sparsity  of  tabulated  values  of  the  functions  (especially  for  large 
values  of  (3  A  ,  as  well  as  the  lack  of  a  thorough  study  of  the  properties 
of  the  functions),  and  the  second  is  the  difficulty  of  handling  the  large 
systems  of  equations  which  arise  in  the  coated  spheroid  problem.  Both 
of  these  obstacles  can  be  surmounted  with  the  aid  of  a  large  scale  digital 
computer.  The  next  few  pages  describe  briefly  the  system  by  which  the 
spheroidal  functions,   V^   and  U^  were  calculated  in  connection  with 
the  present  work. 

The  similarity  of  the  differential  equation,  (37),  for  the  V  -functions 

to  the  corresponding  equation  in  spherical  coordinates,  (12),  suggests  that 
the   V   -functions  can  be  determined  from  a  series  of  the  spherical 
functions,  Thus,  since  the  V5   functions  which  are  the  solutions  to  (12) 
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form  a  complete  orthogonal  set,  it  is  possible  to  represent  the  spheroidal 
functions  in  a  series  as  follows : 

CO 


\=Y,  **~^*  • 


w-1 

If  this  series   is   substituted  into  Eq.    37,  and  the  necessary  conditions  of 

convergence  assumed,   the  result  is 

TV.  "W. 

Equation  12,  with  Jes  —  yiCyi+1),  can  be  multiplied  by  constants  c\fey^ 
and  the  results  added  together  to  give  the  equation 

From  these  latter  two  equations,  the  following  result  can  be  obtained: 

V 
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In  view  of  the  identity, 


V^  =  (1-v^P1 


the  result  in  terms  of  the  more  familiar  Legendre  polynomials  is  equally 
valid : 
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CJL-nC-ft+D-cpJvfjP     =0     •  (49) 


From  this  equation  (for  a  particular  Jk     ),  a  whole  set  of  equations  can  be 
generated  by  utilizing  the  orthogonality  properties  of  the  Legendre  poly- 
nomials, and  this  is  a  particularly  fruitful  step  since,  as  discussed  in 
Appendix  A,  it  results  in  stationary  expressions  for  the  determination  of 
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the  eigenvalues,  as  well  as  recurrence  relationships  for  the  coefficients. 
To  generate  the  new  set  of  equations,  we  can  imagine  multiplying  Eq.  ^9 
first  by  a  Legendre  polynomial  of  order  one,  then  "by  one  of  order  two,  and 
so  on.  Each  equation  so  obtained  is  then  integrated  over  the  complete 
range  of  orthogonality  (  - 1   to  -+- 1  ) .   The  process  can  be  symbolized  as 
follows : 

or 


(50) 


With  the  recurrence  relationship  for  Legendre  polynomials  to  be  found  in 
Section  2.1.3*  "the  last  integral  can  be  transformed  and  evaluated  in  general, 
and  the  value  of  the  first  integral  is  well  known.  The  actual  values  for 
the  integrals  depend  on  whether  the  Legendre  polynomials  are  normalized. 
The  results  for  both  situations  are  included  here,  since  the  coefficients 

o-krx   are  best  obtained  from  the  unnormalized  results,  while  the  eigen- 
values are  most  easily  obtained  from  the  normalized  result. 

If  the  Legendre  polynomials  are  not  normalized,  the  integrals  are 
evaluated  as  follows : 


v/.,    ^   m  Q.Yl+1 
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It  is  clear  from  this  form  that  the  latter  integral  in  Eq.  50  is  zero 
unless  a)  m  =  rv       b ) m=  Yi- a_  or  c)yyv-n+i-  •   Thus,  the  typical 
member  of  the  set  of  equations  has  the  form 
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=  0 


This  is  the  result  which  was  stated  earlier  in  Eq.  39   These  relation- 
ships can  be  employed  to  determine  the  coefficients,  once  the  eigen- 
values, J^k  }   have  been  determined. 

In  the  eigenvalue  calculation,  for  practical  reasons  which  will 
become  clear  below,  it  is  helpful  to  employ  normalized  Legendre  poly- 
nomials.  The  integrals  are  evaluated  in  the  same  way  as  that  indicated 
above  for  the  unnormalized  functions.   The  equations  which  are  so  ob- 
tained can  be  symbolized  as  follows : 
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or  by  the  matrix  equation 
in  which 
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The  conditions  for  the  consistency  of  the  set  of  Eqs.  51  is  the  vanishing 
of  the  determinant  of  the  coefficients  as  follows : 

=  o 

(52) 

The  eigenvalues,  Jkj^       ,  are  those  numbers  for  which  Eq.  52  is  satisfied. 
The  problem  of  eigenvalue  determination  is  thus  in  this  case  equivalent 
to  the  problem  of  finding  the  eigenvalues  of  a  symmetric  (infinite) 
matrix. 

The  main  object  in  developing  the  spheroidal  functions  in  the  fashion 
described  is  to  take  advantage  of  the  service  routines  in  the  ILLIAC 
library  which  are  capable  of  finding  the  eigenvalues  of  symmetric 
matrices  as  large  as  128  x  128,  with  accuracy  of  10  or  11  decimal  places. 
The  fact  (as  is  shown  in  Appendix  A)  that  the  foregoing  relationships 
have  stationary  properties,  together  with  the  potential  sizes  of 
matrices  which  can  be  handled,  makes  it  possible  to  determine  the  eigen- 
values of  nearly  all  interesting  orders  of  the  functions  in  a  re- 
latively short  time. 

In  the  actual  program,  the  equations  (and  therefore  the  eigenvalues) 
were  scaled  by  ( (6J0  ),  and  a  routine  was  written  which  generates  the 
matrix  elements  in  the  form  required  by  the  matrix  eigenvalue  service 
routine.   It  is  worthy  of  special  mention  that  the  equations  (coefficients) 
having   TO.   odd  can  be  handled  independently  from  those  having  7\.  even, 
which  means,  for  example,  that  to  obtain  a  set  of  X   eigenvalues,  one 
needs  to  handle  only  /^     order  matrices.  The  eigenvalues  for  (B>!  —  5,8, 
and  12  are  presented  in  Table  7-  The  effect  of  matrix  size  was  also  in- 
vestigated, with  the  surprising  result  that  for  £1=  8  ,  for  example, 


1+7 
except  for  the  last  five  or  six  values,  a  23  x  23  matrix  gives  results 
which  agree  with  those  given  by  a  to  x  to  matrix. 

With  the  eigenvalues  so  determined ,  Eq  39  can  be  employed  to  find 
the  coefficients.  However,  as  with  the  Bessel  functions  of  the  first  kind, 
(see  Section  2.1.3),  a  straight  upward  use  of  the  recursion  form  is  not 
generally  practicable.  In  the  programs  for  this  work,  the  coefficients 
for  the  lower  order  functions  were  evaluated  by  a  downward  continued 
fractions  technique,  the  form  of  which  can  be  obtained  by  rewriting  Eq. 
39  in  a  fashion  similar  to  that  described  in  Section  2.1.3  in  connection 

A. 

with  the  u-^  functions.  For  the  higher  order  functions,  this  scheme 
was  employed  for  those  coefficients  <5-j^  in  the  range  "A.  =  JfL      .  Then 
a  switch  was  made  and  the  remaining  coefficients  A^^  «>  1  —  Vi  —  Jk.    ,  were 
calculated  from  the  recurrence  formulas  in  a  straight  upward  fashion 
(this  double-barreled  approach  proved  to  be  absolutely  essential  in  the 
calculation  of  the.  higher  order  functions).  The  coefficients  were  all 
generated  on  the  assumption  that  &k!k~   ^"  >   an<i  then  were  normalized  by 
dividing  each  by  dividing  each  by  the  quantity,  ^Aj^,  3-Tx(-y-+  ^  .  This 
normalization  is  equivalent  to  an  adjustment  of  the  functions  so  as  to 
give  the  following  equations 

This  result  is  particularly  convenient  in  the  numerical  work  associated  with 
Eqs.  te  and  tea. 

The  quantities  which  are  needed  in  the  radiation  problem  are  those  as 
follows  (see  Eq.  to  for  example): 
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FIGURE  10   SIMPLIFIED  FLOW  CHART  FOR  MACHINE  CALCULATION  OF  SPHEROIDAL 
V  -FUNCTIONS 


1*9 


These  quantities  are  computed  "by  the  ILLIAC  using  a  program,  a  simplified 
flow  chart  for  which  is  given  in  Figure  10.   The  details  of  the  calcula- 
tion, and  some  of  the  values  computed  will  "be  published  in  a  separate 
technical  report. 

The  next  problem  is  the  evaluation  of  the  functions   U   ,  which 
satisfy  Eq.  36.   It  is  necessary  to  obtain  two  independent  solutions  for 

A  A. 

Eq.  36,  those  which  are  respectively  analogous  to  the  J^  and   iM^ 
functions  in  spherical  coordinates,  together  with  their  derivatives. 
These  independent  solutions  are  identified  by  the  symbols  (Ajg  ±_  and  CAj^ 
in  the  preceding  sections.  The  calculation  of  the  function   lljka.        and 
its  derivative  presents  the  most  difficulty. 

Equation  36  bears  a  strong  resemblance  to  the  radial  differential 
equation  in  spherical  coordinates  (Eq.  11  with  As—  ^Uth-1)  ).  It  seems 
appropriate  therefore  to  try  to  express  the  functions,  \Jj^      ,  in  terms  of 

A 

the  functions,      2         •      Consequently,   it  is  assumed  that  the  solutions 
to  Eq.   36  can  be  expressed  as  follows: 

U,=  L^J^P^    •  (53) 

—  <x> 

Q 

(Publications  by  Schmid  imply  that  such  a  representation  may  be  mathe  - 
matically  proper.)  Assuming  the  necessary  conditions  of  convergence,  this 
series  is  substituted  into  Eq.  36.  The  properties  of  the  function  as  dictated 
by  Eq.  11  are  employed  to  reduce  the  result  to  the  statement 
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With  a  repeated  use  of  the  expression  for  the  derivative  function, 

the  term  involving  the  second  derivative  can  he  replaced,  and  Eq. 
5  4  put  into  the  form 

OO  A  A 

<=ao  A, 

A 

If  now  the  coefficients  of  the  terms,  J^ ,  are  individually  put  equal 
to  zero,  the  resulting  expressions  can  he  put  into  a  form  as  follows: 


Ca.-n.-3X  sen- 1) 


H- 


(55) 


Note  the  similarity  of  this  relationship  to  that  of  Eq.  59-  The  simi- 
larity is  such  that  if  the  condition  CL^  —  <&fejt  were  imposed,  the  magni- 
tudes of  the  remaining  coefficients  would  he  equal.   (The  signs  would  he 
such  that  <Xh  ^   ^=1  — dljj  ^    ,  G^  .  =  A*  ^+4  >   and  so  on>  alternating.)  As 
with  the  coefficients,  if  Jk  is  odd  (even),  only  the  odd  (even)  co- 
efficients are  non-zero. 
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A  glance  at  the  first  term  in  (55)  shows  that  the  positive  side  of 
the  recurrence  chain  breaks  off  at  Y\.=  l   and  >x=  a,   ,  while  from  the 
last  term  it  is  clear  that  the  negative  side  of  the  recurrence  chain 
"breaks  off  at  Y\.  =  -a.   and  ~a=-3   .  As  a  result,  the  value  of  CL^ 
corresponding  to  one  of  these  particular  four  values  can  be  assigned, 
and  on  the  basis  of  this  assignment,  a  positive  (negative)  side,  odd 
(even)  set  of  coefficients  can  be  determined.  Furthermore,  if  it 
happens  that  Q-j^i.  =  ^M.-^  *  °^,  a,—  Q-js  -3  >   "then  the  coefficients 
generated  on  the  positive  side  from  Cl^  ^_   (  o~^^  )  are  the  same  as  those 
generated  on  the  negative  side  by  CLj^  (  <Xj^         )}   since  in  Eq.  55,  tl 

can  be  changed  to-trrL+1}  without  a  change  in  the  equation.  Thus,  suppose 
for  example  one  of  the  1  coefficients,  say  <Xkl     ,  is  assigned  a  value. 
Then,  since  Eq.  36  is  a  differential  equation  of  second  order,  there  is 
still  one  constant  to  be  assigned*   If  this  second  constant  is  set  by 
putting  (X.   — O  C  Q-Jk,\.      assigned),  the  coefficients  are  all  determined, 
and  the  series  gives  the  spheroidal  functions  of  the  first  kind,  as 
follows : 

-*.=  -1,0 

where  the  prime  indicates  that  the  sum  is  on  odd  (even)  values  of  tl   only. 

(The  even  set  results  from  the  assignment  of  ^>+^  >   with  (Xj^  _.—  o   .) 
On  the  other  hand,  if  <X^  _x  (Ol^.^^  )  is  assigned  to  fix  one  of  the  con- 
stants, while  a*ia  (  CC*  1  )  is  put  equal  to  zero  to  fix  the  other,  the 
series  is  /<rj\ 


\,=YL'  ^A^  • 
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In  practice  it  is  convenient  to  take  advantage  of  the  relationship 

a  A 

j,  =  c-ir  n 

--IV.  rv-i_ 

and  so  transform  the  series,  (57)  >  into  a  sum  over  positive  values  of 
TV  , 

^-ip  (58) 

in  which  the  coefficients  are  identical  to  those  in  Eq.  56   Unfortunately, 
the  series  (58)  converges  rather  slowly  when  7a.   is  in  the  neighborhood 
of  unity.  Some  difficulties  were  encountered  in  the  programs  which  were 
written  for  the  ILLIAC  in  the  calculations  of  the  higher  order  functions, 

A 

since  the  magnitudes  of  the  P"v  -functions  approach  machine  capacity 
before  the  functions  converge  to  a  desirable  accuracy. 

A  simplified  flow  chart  of  the  program  to  calculate  the 
-functions  is  given  in  Fig.  11.  An  interesting  feature  of  this  pro- 
gram is  the  print  out  which  indicates  that  there  is  an  insufficient 
number  of  terms  available  and  also  provides  a  measure  of  the  over -all 
accuracy  of  the  calculation.  This  latter  is  an  internal  check  which  is 
based  upon  the  value  of  the  Wronskian  determinant  of   Uj^i.  and  Uj^a.  • 
This  quantity,   \A/  t   has  a  constant  value  as  can  be  shown  as 
follows : 

<a.u. 


5U 

Multiply  the  first  of  these  latter  two  equations  by   ^Jfe,^  ,   the  second 
by   CAj^  ,  and  subtract,  to  obtain  the  result 

or 

~\y\7~   -=.      constant   . 

The  value  of  the  constant  can  be  found  by  obtaining  the  solutions  to 
(36)  in  the  limit  as  "U_   approaches  infinity.   In  that  case 

U  =  A  coS(3JHl  ■+  B  SlTfl  Piu 
so  that 

W=  -pi  • 

In  the  program,  after  the  independent  calculation  of  each  of  the 
functions   (i^  ,      LLkL    ,    u^  ,  and   tXfc,,  ,   the  quantity 

is  formed  in  the  machine  and  the  result  printed  out  along  with  the  functions. 
The  closeness  of  the  number  printed  to  the  exact  value,  -1  ,  is  a  good 
measure  of  the  over-all  accuracy  of  the  calculation.  A  facsimile  of  a 
typical  print  out  is  given  in  Table  2. 
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TABLE  2.   FACSIMILE  OF  A  TYPICAL  ILLIAC  PRINT  OUT 
OF  THE  SPHEROIDAL  FUNCTIONS,   L) 


The  first  two  lines  of  print  indicate  that  there  were  not  enough 

A  A  . 

of  the  functions   N^  ,  and  (ni^  to  calculate   U33L   and  il3-L      to 
the  specified  accuracy.   In  these  first  two  lines,  the  first  number  is 
the  value  of  the  last  term  computed,  while  the  second  number  is  the  cumu- 
lative sum  (both  scaled).  These  numbers  indicate  that  in  this  case  the 


functions 


U 


and 


u: 


3a,  are  good  to  about  seven  figures,  and  this  is 


born  out  by  the  amount  by  which  the  scaled  Wronskian  differs  from  -1 


Number 


exp. 


Number 


exp. 


+21992 

-281529 

=69 
-68 

+2061J+ 
+127821 

-62 
-61 

u 

W/f3£ 

+11000 

+01 

-10000071 

+01 

u 

X. 

i 

u 


•u' 


+7^9599554   +oo  -162^53319  +01  +206136057  +00  +10225702^  +02 
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For  convenience  in  the  calculation  of  the  radiation  pattern,  the  arbi- 
trary constant  (  Ctj^j.  and  so  on)  can  be  fixed  in  such  a  way  that  as  PH.- — >°° 

f 

To  accomplish  this  normalization,  it  is  only  necessary  to  enforce  the 
following  condition: 

or 

CO 

E  ^       =  1  • 

Since  the  sum  is  over  only  odd  (even)  tl  ,  this  is  accomplished  in  the 
machine  by  making  O^fe^  positive,  and  changing  the  signs  of  alternate  co- 
efficients in  the  summation  process. 

The  radiation  patterns  of  the  uncoated  spheroid  are  then  easily 
calculated,  with  the  spheroidal  functions  determined  in  the  way  described 
above.  The  pattern  calculation  is  based  on  Eqs.  k2   and  k^>   and  is  done 
almost  exactly  as  described  in  Section  2.1-3,  with  the  spheroidal  functions 

A  ,  (  A  , 

replacing  the  spherical  functions  as  follows:    ilj^t  *~~*  L-'n,  ■>  ^Aa.*  "'  ^vo 
and    y^tl-V*}      <^-^    \\^        .     Some  of  the  results  are  presented  in 
Figs.  21  to  27,  with  some  corresponding  coated  patterns  alongside  for  com- 
parison. The  patterns  in  Fig.  21,  for  the  smaller  values  of  ^   ,  were 

calculated  on  a  desk  calculator  using  functions  modified  from  those  computed 

q  10 

by  Myers  and  Wells  ,  and  coefficients  tabulated  by  Stratton,  et  al.     These 

latter  calculations  were  done  early  in  the  investigation  to  study  the 

feasibility  of  the  method  and  are  subject  to  greater  error  than  the  others; 
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in  these  calculations,  the  digital  computer  was  employed  only  to  solve  the 
systems  of  equations  which  arise  in  the  coated  spheroid  problem  (12 
equations  in  12  unknowns  was  the  maximum  size  in  this  case). 

The  biggest  job  in  the  computation  program  is  in  connection  with 
the  coated  spheroidal  radiator.  The  calculation  is  based  on  Eqs.  k'J   and 
h8   and  the  operations  indicated  therein  are  carried  out  in  exactly  the 
manner  indicated.  Fortunately,  the  odd  and  even  orders  can  be  handled 
separately.  The  pair  of  equations,  (V7)  and  (hS) ,  are  replaced  by  the  four 
equations  obtained  by  equating  separately  the  real  and  imaginary  parts  of 
(hf)   and  (k&) .     A  guess  is  made  as  to  the  number  of  equations  required  for 
convergence,  and  the  augmented  matrix  for  this  number  of  equations  is 
generated  by  the  machine  in  the  form  required  by  the  ILLIAC  library  service 
routine  which  obtains  the  solutions  of  a  set  of  linear  equations  (a  max- 
imum of  1^3). 

The  variation  of  the  coefficients,  C;  ,  as  a  function  of  the  number  of 
equations  was  investigated.  The  results  of  one  such  investigation  are  pre- 
sented in  Table  3.  It  is  clear  from  these  data  that  the  process,  as  dis- 
cussed near  the  end  of  section  2.2.2,  is  actually  a  convergent  one,  and  in 
principle  could  be  carried  out  to  any  desired  accuracy.  The  coefficients, 

Re.  C,    ,   and  "X*a  C,  ,    so  calculated  are  fed  back  into  the  machine  and 
a  3 

used  in  the  basic  routine  for  radiation  pattern  calculation. 

We  now  leave  the  study  of  prolate  spheroidal  radiators,  and  take  up 
once  more  the  study  of  the  sphere,  in  order  to  remove  the  restriction  of 
axial  symmetry  on  the  sources  . 
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TABLE  3.   COMPARISON  OF  THE  COEFFICIENTS  CALCULATED  FROM  DIFFERENT  NUMBERS 
OF  EQUATIONS  IN  THE  COATED  PROLATE  SPHEROID  PROBLEM. 


The  similarity  in  the  numbers  shows  that  the  infinite  system  can  he  re- 
placed by  a  finite  system  of  equations.  The  example  is  for  the  parameters 
as  follows:    $,1=5",  (3.JL  =  8  ,  ux  =  l.OTT  ,  7^=  l.i  ^  slot  oJT  V  =  0. 


20  equations 

24  equations 

R-S 

iQ   equations 

32  equations 

1 
5 
5 
7 
9 

.035669646 
-.061748146 

.009879699 
-.000092924 

.000000405 

•  035669590 
- . 061748740 

.009881601 
-.000092948 

.oooooo4o6 

.035669589 
-.061748754 

.OO9881658 
-.000092949 

.000000406 

.035669589 
-.06174^750 

.009881652 
- . 000092947 

.000000405 

11 

-.000000001 

-.000000001 

- . 000000001 

13 

. 000000000 

.000000000 

15 

ImCj 

.000000000 

1 

3 

5 

7 

9 

11 

13 

15 

.048023926 
-.009061661 

.033134300 
-.002072803 

.OOOO76310 

.048023784 
- . 009061355 

. 033143465 
-.002083607 

.000083883 
-.000002187 

. 048023780 
-.0Q906l3^5< 

.033143727 
- . 002083989 

. 000084252 

- . 000002383 

.000000045 

. 048023780 
- . 009061346 

.033143723 
- . 002083994 

.000084262 
- . 000002392 

.000000048 

- . 000000001 
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3.   SPHERE  WITH  ARBITRARY  SOURCE  CONFIGURATIONS 

The  analysis  for  the  fields  which  result  from  the  excitation  of  a 
spheroidal  conductor  by  an  arbitrary  slot  is  very  complicated.,   In  fact, 
it  appears  that  only  the  sphere  can  be  treated  numerically  with  existing 
techniques.   However,  in  spherical  coordinates,  it  is  possible  to 
represent  the  fields  in  terms  of  a  pair  of  scalars,  as  indicated  in  the 
following  developments.  In  this  process,  one  or  the  other  of  the  field 
variables  is  eliminated  from  Maxwell's  equations  to  obtain  the  vector 
Helmholtz  equation 

VxV*  B  -  u/Ve  B  =°         (59) 

where  Q  is  either  E  or  H   «  The  problem  is  then  simplified  by 
the  division  of  all  possible  fields  into  parts  in  which  either  E  or 

H  is  transverse  to  the  radial  direction,  with  consequent  representation 
in  terms  of  a  vector  (or  Hertz)  potential  which  has  only  an  -f  -component. 
Thus,  for  TM  fields 

H  =  y*f+-  (6o) 

and  for  TE  fields 

E=V*3^  (6i) 

The  subtitution  of  (6o)  into  (59)  gives  for  example 

yxVxfr-^t  =  vV  (62) 

where    V     is  arbitrary.      If    ^     is   chosen  so  that     W  =  JiL  the 

following  equation  is  obtained  for       r         : 


6o 


(63) 


and  a  similar  result  can  be  obtained  for  Q    .  The  scalar  Helmholtz 
equation,  (62),  is  separable  in  spherical  coordinates,  and  the  solutions 
have  the  form 


4^H  =-^v.R"vT"  +^'vT" 


(6k) 


with  the  constants  selected  so  that  the  boundary  conditions  are 
satisfied.  For  the  moment,  for  brevity,  let  the  solutions  for 
be  represented  as 

f  =  (&rtT'ce,+)  (65) 

where  \\Q*)   and  I  (9>4>}  are  defined  appropriately.  With  this 
notation,  expressions  for  the  fields  can  be  written  as  follows : 

a)  TM  fields 

^&E  =  -^-R'vT  +^'vT'  (6?) 

b)  TE  fields      C3  =  R.V>  T  "(.  9^ 

£   =  v*  R"T"r  ■=  -t^'vT"  (68) 


(69) 
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Thus  for  example,  the  transverse  component  of  any  electric  field  can 
be  written 

ET  =  YL  R!i,  vT,"  *  £  +  aft  yt  (to) 

where  the  constants  in   R.   ,    \    ,   f^"   ,  and   ry1'   ,  are 
determined  by  application  of  the  boundary  conditions,  and  the  summation 
over  the  index   i-   is  to  be  interpreted  as  a  summation  over  terms 
involving  all  permissible  combinations  of  V   and  ~rr^     in  Eq.  64. 

In  any  complete  region  (including  both   B  —  O   and  0  —  Tf   ); 
the  constant   cXy^     in  Eq.  64  must  be  zero,  in  order  to  avoid 
infinities  on  the  z-axis,  and  the  index   /    must  be  an  integer.  To 
simplify  the  analysis,  it  is  convenient  to  treat  separately  the  parts  of 
the  solution  involving  snrbTm,^   (odd  solutions)  from  those  involving 

cos  -m.4>      (even  solutions).  It  can  also  be  seen  that  the 
quantity  \\_      ,  as  defined  from  Eq.  64  will  contain   S~    as  a 
multiplier.  It  is  convenient  therefore  to  represent  the  radial  vari- 

A  A 

ations  in  terms  of  the  functions   J^)  '—  **  V  &>         and  Nl^  =  ^"V^ 
which  were  employed  in  the  analysis  in  the  previous  sections  (let  the 
constants  absorb  the  quantity   p   which  is  also  needed  as  a  multiplier  in 
order  to  change  the  functions).  With  these  changes,  and  the  definitions, 

TLt   =  PC  cos  -mc{>    and  T10  —  f^  >  Sm  yti$  ,  Eqs .  70  take  the 

following  form: 

+  c  <_c  [L  +  K^,  N^  VX0  *  >r  (71) 

+  jnkL<™  aJ-  +  A—*  ^N^  vTU 

3  &T"  "XT 

4   fc  cJL^r  cJL-T 
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It  is  convenient  in  the  theoretical  discussions  to  define  vector 
functions  as  follows : 

el?  =  S?X%  ^  r     ^   T/s  =  vXc 

along  with  appropriate  scalar  functions  of   -f-    so  that  Eq.  71  can  be 
written  in  the  simpler  form 

ET  =  XL  OP^*-  +%:^eZ  4-^Cf^  +t(^S   ,    (72) 

It  is  shown  in  Appendix  B    that  the  vector  functions  e^e 
and    e«e     have  orthogonality  properties  which  permit  field  expansions 
of  the  vector   Er   rather  than  one  component  of  it  as  was  done  in  the 
case  of  axially  symmetric  excitation.  In  the  next  section,  the  use  of 
these  expansions  in  the  solutions  for  a  conducting  sphere  in  free  space 
is  illustrated.  The  development  serves  the  dual  purpose  of  providing  a 
review  for  the  reader  and  standardizing  the  notation  so  that  comparisons 
to  the  dielectric  coated  sphere  problem  can  be  more  readily  made0 

$.1  The  Slot  Excited  Conducting  Sphere  in  Free  Space 

Let  a  conducting  sphere  of  radius  -^   be  excited  by  applying  a 
(known)  electric  field  in  an  opening  which  is  cut  somewhere  on  its 
surface.  Then,  the  electric  field  is  known  everywhere  on  the  surface, 

-t*  =  -fx         ,  viz.,  it  is  zero  over  most  of  the  surface  and  ET<X, 
in  the  slot  opening.  Thus,  from  (72),  the  following  equation  can  be 
written: 

E_  I   =-E.  -  7L  XU&^eSi  -4-  X:c^)e!i  +%l(fc)e^  +  %'tf^K  Wi 
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Since   ETo,  is  known,  the  orthogonality  relations  listed  (and 
demonstrated)  in  Appendix  B    can  be  employed  to  determine  the 
following  equation  for  the  excitation  expansion  coefficients: 

TCcpwirt   =    Lr  E^  el  AS  (7^) 

Tk'cfStirt   =  Lr  ETe;  e,:  AS  (75) 

X=y-.  e£a£S 

Ulcp^O  »  X,o7  EtJ  e4  AS  (76) 

L.  e£  &s 

TOfW  =  LrE^-e^AS  (77) 

Of  course,  Eq.  73  with  the  coefficients  calculated  from  (7*0- 
(77)  must  he  exactly  equal  to  Eq.  71  with  -T   set  equal  to   -rj 
For  the  problem  of  the  conducting  sphere  in  free  space,  however, 
Eq.  71  is  simplified  because  of  the  requirement  that  the  solution  must 
represent  outgoing  waves  as  -^  — >  00    .  This  implies  that  the 
constants  in  (71)  must  be  related  as  follows: 

0-™.  =  ~i  Kir.       ,    Cn^  =  ~ 4  ^  -  (78) 

Thus,  the  expression  for  the  transverse  electric  field  can  be  written 


6k 


E T  =  YL  <*4  H^p.^  e.'s  -h  c^  dL  HACM  e/s    .     (79) 

*"  ^£»  At 

An  evaluation  of  this  expression  at  -f-  =  -ft      ,   and  subsequent  comparison 
of  it  to  Eq.  73  indicates  the  manner  in  which  the  constants  can  be 
evaluated  numerically,  and  the  result  is 


XUt.)     =  Lr  e^  e4  <*-$ 


(80) 


U  =  rf      -  *gj=l£flig  •  (81) 


— — — ^—  u-fz.T 

JL-r  A* 

The  best  form  for  the  computation  of  the  fields  is  probably 

«. — r  A  f  i;         ^  r  cos  w4>"i      £         r-^t^r-sivt^^T 


jjtu&o        ^-T 


^  -^q-  -|-  sm 9 


since 


and 


c^  *  "Xe" 


-tsm-t 


A  ^-m.  C—  Siw.™.*  -)  A        n^  f  COS  Tn(j)  ") 

r;  =  7T«  ^  =  1_™_  RTi  c*s  m*  J  -4.  &?»  1  ^m*  5 
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Also,  in  view  of  these  latter  two  relations,  the  following  forms  for 
the  evaluation  of  the  integrals  in  (7^)  -  (77)  can  be  written  down: 
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FIGURE  12  SLOT  EXCITED  SPHERICAL  RADIATOR  WITH  DIELECTRIC  COATING 
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3.2  The  Slot  Excited  Conducting  Sphere  Coated  with  Dielectric 

When  the  conducting  spherical  antenna  is  coated  with  dielectric 
material,  the  exterior  space  is  split  into  two  homogeneous  regions,  the 
region  (i)  of  the  spherical  shell  of  dielectric  (  +,  —  ""f"  —  ^\  ) 
and  the  region  (il)  exterior  to  the  shell  (    -f  >  -T^    ).  As  in  the 
problem  of  the  conducting  sphere  in  free  space,  the  tangential  electric 
field  is  presumed  to  "be  zero  on  the  surface   -f  —  -tj   ,  except  over 
the  slot,  where  it  has  the  prescribed  value    ET(v    .In  the  region 
(i),  the  transverse  electric  field  has  the  form  precisely  as  given  in 
(7l)«  The  field  over  the  conducting  surface,  -T  —  ^  ,  can  be 
represented  in  the  form  of  Eq.  73*  in  which  the   If      's  are  regarded 
as  the  excitation  coefficients,  determined  from  relations  (7*0  -  (77) « 
The  expression  given  in  (71)/  evaluated  at  -f  ■=•  "Tj1  ,  must  be  equal 
term  by  term  to  an  expression  of  the  form  of  (73)/  to  satisfy  the  boundary 
condition  at  the  conducting  sphere.  This  implies  the  following  set  of 
four  equations: 

/v  a 

a-L.   XC&>    ■+    -CSN.A^  =  Xl%  (82) 


(83) 


Next,  the  boundary  conditions  at  the  outer  surface  of  the  dielectric  shell; 

-T  =  -^  are  applied.   In  region  (il),  "t*  >  *t^  ,    the  transverse 
electric  field  has  the  same  form  as  that  given  in  (79)/  since  it  must 
consist  of  outgoing  spherical  waves  as  ^  >co    ,  When  the  condition 
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ETr    —  L  Tir      is  enforced  (where  again  the  form  of  Eq.  71 
for  the  field  in  region  (i)  is  employed),  the  result  is  the  equation 


LUu.-i^  +C„N,(p. 


-  // 


1  4 


— "A+        (85) 


This  equation  is  satisfied  if  the  following  set  of  four  equations  is 
satisfied: 

<^.Jjfito  +  *lA&*  =  i"*Rj|^  (86) 


"  ~£^~        °  A+-       £°  °    A-t- 

The  boundary  condition  on  the  magnetic  field   H  _1   —  H. 


(87) 


must  now  be  enforced.  A  careful  examination  of  the  equations  (66)   -   (69) 
shows  that  an  equation  for     HTz    can  te  bitten  as  follows : 

,.  1  u;:e  &jcw  +  41.  an, 


(88) 

A  A 
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In  region  (il),  the  functions  must  represent  outgoing  traveling  waves. 
Thus,  the  boundary  condition  on  the  magnetic  field  gives  the  following 
equation: 

»     i»>*>-  TT  JL+  (89) 

=  7  -J,  {IsiH^^vT.  +  ^H^WT^*  . 

(The  terms  have  orthogonality  as  mentioned  under  Eq.  72,  where  the 
notation  was  used  for  the  vector  functions.)  This  equation  is  satisfied 
if  the  following  set  of  four  equations  is  satisfied : 

<„s  AjX<p.«  +  jC*aN^  =  **  tl»4K*£wa    (90) 


The  forgoing  equations  imply  that  the  TM  and  TE  field  coefficients  can 
be  handled  separately.  The  TM  coefficients  can  be  found  from  the  sets  of 
Eqs.  83,  87,  and  SI.  These  sets  of  equations  are  regarded  as  independent 
sets  of  equations  for  the  field  coefficients,  and  their  solution  is  indi- 
cated by  the  display  of  the  following  matrix  (the  augmented  matrix  of  any 
particular  set) ; 
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The  TE  coefficients  can  be  found  from  the  sets  of  EqSo  82,  86, 
and  90.  The  augmented  matrix  of  any  one  set  of  these  coefficients  is 
as  follows : 


O 

A 


A- 


cL^- 


Mr 


dL^ 


if 


O 


TE- 


For  the  radiation  pattern  calculation,  the  coefficients  JL 
and  q  are  needed  explicitly.  The  results,  obtained  frc 
foregoing  equations,  can  be  written  compactly  as  follows: 


o   °      u  2 

ATM)t 
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(92) 


nTn_e 


A. 


(95) 


where   ZATf^   and   Zate   are  respectively  the  determinants  of  the 
first  three  columns  of  the  TM  and  TE  matrices  which  are  written  down  above. 
Note  that  the  determinant  for  the  TM  coefficients  is  identical  to  the 
corresponding  determinant  (Eq.  29)  in  the  case  of  the  symmetrically 
excited  sphere.  Hence,  the  numerical  results  presented  earlier  can  be 
applied  here. 

As  in  the  case  of  axially  symmetric  excitation,  the  effect  of  the 
dielectric  coating  can  be  assessed  by  examining  the  ratios  of  the 
corresponding  expansion  coefficients  in  the  coated  and  uncoated  problems : 

3^8   =  £r  ft    AhyfrO  (<*) 


to^        AT6„  (95) 


The  ratio  in  Eq.  9^  was  studied  in  Section  2.1.2;  numerical  results  are 
given  in  Figs.  ^+0  to  63,  and  Table  6,  The  ratio  in  Eq.  95  has  not  been 
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studied  numerically,  but  in  view  of  the  fact  that  the  only  difference 
from  the  former  is  in  the  replacement  of  some  of  the  functions  by  their 
derivatives  (and  the  relationship   £_£H  =  -£±±   Z^  —  ZL  +  A    ), 

OLTC         a-  — 

it  is  likely  that  the  latter  ratio  exhibits  a  similar  behavior. 
Unfortunately  the  numerical  study  of  the  multitude  of  problems  which 
are  represented  formally  by  the  analysis  (given  in  the  preceding  pages) 
of  the  coated  spherical  radiator  with  arbitrary  slot  configurations 
must  be  postponed  to  later  investigations. 
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h.      SUMMARY  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

Analytical  methods  have  been  presented  for  the  determination  of  the 
changes  in  the  radiation  patterns  brought  about  by  dielectric  coatings  on 
spheroidal  radiators.   Numerical  results  are  included  for  spheres  and 
prolate  spheroids  which  are  excited  by  axially  symmetric  slots .  A 
description  is  given  of  the  machine  computation  system  which  makes  the 
analytical  results  useful.  This  system  includes  a  program  for  the 
semi-automatic  calculation  of  the  radiation  patterns  of  spherical 
radiators,  a  practical  method  for  the  determination  of  the  eigenvalues 
of  all  of  the  orders  of  the  spheroidal  functions  which  are  needed  in  order 
to  obtain  convergent  representations,  and  a  new  method  for  the  calculation 
of  the  spheroidal  "radial"  functions. 

The  comparison  of  the  radiation  patterns  of  an  uncoated  spherical 
radiator  to  those  of  a  dielectric  coated  radiator  having  the  same 
excitation  shows  that  the  effects  of  the  coating  may  be  very  slight  in 
some  cases  and  very  important  in  other  cases.  One  effect  of  a  coating  of 
the  order  of  several  hundredth's  of  a  wavelength  on  a  moderate  size  radiator 
(a  few  wavelengths  in  diameter) is  to  introduce  into  the  patterns  some  qf  the 
characteristics  of  a  radiator  which  is  somewhat  larger,  with  an  increase  in 
the  radiation  in  the  general  direction  of  the  poles.  In  the  patterns  in 
which  the  excitation  is  by  means  of  a  nonequatorial  slot  (Figures  15  and 
16,  for  example),  the  increase  in  radiation  in  the  general  direction  of  the 
opposite  pole  is  similar  to  that  which  might  be  expected  in  the  event  of  an 
increase  in  surface  type  waves  traveling  around  the  sphere,  but  the  theoretical 
representations  do  not  obviously  permit  such  an  interpretation.   On  the  other 
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hand,  as  indicated  in  Figures  17  to  20,  as  the  permittivity  or  coating 
thickness  is  varied,  the  patterns  may  go  through  the  successive  stages 
of  slight  deterioration  from,  complete  deterioration  from,  and/or 
smoothing  of  the  patterns  of  the  uncoated  radiator.  If  the  quantities 
R  and  I   (  pg.  26  )  which  appear  in  the  expressions  for  the  coated 
spherical  radiator  are  studied  in  the  limit  as  the  dielectric  constant 
tends  to  infinity,  it  can  "be  seen  that  the  quantities  in  the  symbolic 
brackets  tend  to  sines  and  cosines,  and  become  independent  of  n.  Thus, 
except  in  the  event  of  a  zero  in  the  argument  of  the  ^cosines;  (f3,  (rp-r  )), 
the  shape  of  the  pattern  tends  toward  that  of  a  conducting  sphere  having 
a  radius  equal  to  the  outside  radius  of  the  dielectric  casing.  This 
effect  is  apparently  beginning  to  show  up  in  the  patterns  of  Figure  18 
for  the  higher  dielectric  constants,  although  the  smoothing  is  somewhat 
surprising.         .   < 

In  retrospect,  it  appears  that  -sameoof. (the  "effects  :might  be 
better  understood  if  the  analysis  were  recast  in  terms  of  spherical 
transmission  line  theory.  In  the  latter  theory,  the  cpntributions 
to  the  total  field  made  by  the  different  orders  of  the  functions  are 

1  regarded  as  modes  propagating  along  a  set  of  independent  transmission 

11 
lines.  That  such  a  representation  is  possible  has  been  shown  by  Marcuwitz 

12 
Felsen,  and  others   .  In  principle  one  could  calculate  a  set  of  mode 

transmission  coefficients  to  relate  the  excitation  coefficients  (or  mode 

input  voltages)  to  the  distant  field  coefficients,  (or  mode  output  voltages) 

;  which  appear  at  the  output  of  a  kind  of  spherical  transmission  mode  filter. 
The  radiation  patterns  of  a  few  prolate  spheroidal  radiators  are 

j  presented.  Unfortunately,  the  system  of  calculation  employed  to  obtain 
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these  patterns  is  very  much  more  clumsy  and  time  consuming  than  the  sphere 
calculations.   However,  the  calculations  which  were  carried  out  in 
connection  with  the  coated  spheroid  problem  are  interesting  for  reasons 
which  extend  beyond  the  specific  problem.  For  it  has  been  demonstrated 
that  the  difficulties  (of  enforcing  the  boundary  conditions)  which  arise 
in  problems  in  which  the  (eigen-)  functions  are  different  on  two  sides  of 
a  boundary  may  not  be  insurmountable.  The  results  in  Table  3  indicate  that 
satisfactorily  accurate  results  can  be  obtained  by  including  only  a 
manageable  number  of  equations  out  of  the  infinite  set  which  is  required 
for  an  exact  representation. 

Comparison  of  the  patterns  of  the  prolate  spheroidal  radiators  to 
those  of wthe  spherical  radiators  lends  support  to  the  idea  that  the  effects 
associated  with  spheroidal  radiators  resemble  those  associated  with 
spherical  radiators.  This  similarity,  coupled  with  the  difficulty  in 
carrying  out  the  prolate  spheroidal  solutions  provides  further  motivation 
for  studies  of  the  spherical  radiator,  in  order  to  be  able  to  better  judge 
which  of  the  spheroid  problems  are  likely  to  be  the  more  interesting. 

The  following  additional  work  is  recommended: 
l)  The  dielectric  coated  sphere  problem  should  be  studied  analytically 
and  numerically  from  the  point  of  view  of  spherical  transmission  line 
theory  and  the  feasibility  of  dielectric  layers  as  spherical  transmission 
line  mode  filters  should  be  investigated.   2)  The  results  of  the  present 
analysis  for  spherical  radiators  with  arbitrary  slot  configurations  should 
be  applied  to  such  problems  as  a  traveling  wave  slot  on  a  sphere.   3)  The 
method  of  calculating  the  spheroidal  functions  should  be  extended  and 
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improved  since  the  method  has  promise  in  many  other  problems  (for 
example,  slot  excited  radiators  whose  cross  section  has  an  abrupt 
change  in  size  and  character) . 
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TABLE  4.   EXCITATION  COEFFICIENTS  FOR  AN  AXIALLY  SYMMETRIC  SLOT  IN  A 

SPHERICAL  RADIATOR 
The  quantity  tabulated  is 


3.TI+1 


r  p^a 


It  is  evaluated  by  a  five  point  Simpson's  rule  integration  over  a  slot 
which  subtends  two  degrees.  Values  presented  here  are   oC_^(  10°)  ^  -yu 
from  one  to  eighty.  Data  and  output  tapes  for  slots  at   9=  \0°(.\O°)R0  • 
Tt :  l(l)80,  are  on  file  at  the  U.  of  I.  Antenna  Laboratory. 


OC. 


l(n)27 

-79192  -03 
+12995  -02 
-17770  -02 
+22137  -02 
-25995  -02 
+29261  -02 
-31864  -02 
+33751  -02 
-34888  -02 
+35260  -02 
-34869  -02 
+33739  -02 
-31912  -02 
+29446  -02 
-264l4  -02 
+22903  -02 
-19010  -02 
+l484l  -02 
-10506  -02 
+6118+  -03 
-17886  -03 
-2375+  -03 
+62729  -03 
-98125  -03 
+12915  -02 
-15513  -02 
+17556  -02 


28(n)5+ 

-19009  -02 
+19855  -02 
-20092  -02 
+19736  -02 
-18817  -02 
+17383  -02 
-15^91  -02 
+13211  -02 
-10623  -02 
+78103  -03 
-48628  -03 
+18702  -03 
+10789  -03 
-38997  -03 
+65138  -03 
-88508  -03 
+10851  -02 
-12465  -02 
+I3658  -02 
-l44o6  -02 
+1^702  -02 
-145^9  -02 
+13967  -02 
-12984  -02 
+11641  -02 
-99907  -03 
+80896  -03 


55(n)80 

-60023  -03 
+37965  -03 
-15418  -03 
-69286  -04 
+284ll  -03 
-48410  -03 
+66^>66  -03 
-81800  -03 
+94320  -03 
-10363  -02 

+10955  -02 
-11199  -02 
+11099  -02 
-IO667  -02 
+99254  -03 
-89065  -03 

+76489  -03 

-61976  -03 

+46023  -03 
-29156  -03 
+11912  -03 
+51733  -04 
-21587  -03 
+36850  -03 
-50534  -03 
+62269  -03 
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TABLE  5.   RESONANCE  AND  CONVERGENCE  FACTORS  FOR  SPHERICAL  RADIATORS. 
The  values  tabulated  are  different  orders  of  the  combinations  of  the  spherical 
Bessel  functions  which  appear  in  the  series  expansions  of  the  fields  of  a 
spherical  radiator  whose  radius  is  2.^  wavelengths.   Data  and  output  tapes  for 
these  quantities  for  different  sphere  sizes  as  follows:  ~t~/?\  —  .l(  .05)1*05, 
tA  =  l.l(.l)2.^V?i=  2.4(.8)4.8,  are  on  file  at  the  U.  of  I.  Antenna 
The  machine  was  instructed  to  calculate  successively  higher 
until  the  exponent  of  the  quantity       N>% 
(Note  +123  -02  means  .00123)  ^  +  ^ 

A.  A 

j:  n 


LaDora 
orders 

cory.   xue  maunxm 
of  the  quantitiei 

-5 

became 

less  than  10 

1 

jAr 

1 

+10044  +01 

2 

+10134  +01 

3 

+10273  +01 

1+ 

+10467  +01 

5 

+10725  +01 

6 

+11062  +01 

7 

+11496  +01 

8 

+12056  +01 

9 

+12790  +01 

10 

+13769  +01 

11 

+1511^  +01 

12 

+17019  +01 

13 

+19758  +01 

11+ 

+23231  +01 

15 

+23902  +01 

16 

+13905  +01 

17 

+40526  +00 

18 

+88558  -01 

19 

+17242  -01 

20 

+30354  -02 

21 

+48o4o  -03 

22 

+68263  -04 

23 

+87336  -05 

24 

+10106  -05 

25 

+10627  -06 

26 

+10202  -07 

27 

+89809  -09 

28 

+72778  -10 

29 

+54488  -11 

30 

+37816  -12 

T 

J'" 


n!7 


+52857  +00 
+91823  +00 
-22277  +00 
-10219  +01 
-39022  +00 
+73515  +00 
+10276  +01 
+29563  +00 
-68837  +00 
-11699  +01 
-96790  +00 
-35622  +00 
+29995  +00 
+75329  +00 
+79978  +00 
+36431  +00 
+69952  -01 

+88559  -02 
+90006  -03 
+75697  -04 

+52945  -05 
+31008  -06 
+15346  -07 
+64788  -09 
+23545  -10 
+74262 
+20477 
+49690 
+10675 

+20412 


-12 
-13 
-15 
-16 

-18 


A. 

A.  , . 

J> 

■N'* 

+85147 

+00 

-41257 

+00 

-98875 

+00 

-48556 

-01 

+95930 

+00 

+75214 

+00 

-30595 

+00 

-10575 

+01 

-89730 

+00 

-90763 

-01 

+75798 

+00 

+12550 

+01 

+13732 

+01 

+13250 

+01 

+13231 

+01 

+11215 

+01 

+63274 

+00 

+29746 

+00 

+13131 

+00 

+55095 

-01 

+21918 

-01 

+82621 

-02 

+29553 

-02 

+10053 

-02 

+32599 

-03 

+10101 

-03 

+29968 

-04 

+85310 

-05 

+23343 

-05 

+61495 

-06 
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TABLE  6.   NUMBERS  WHICH  INDICATE  THE  RELATIVE  IMPORTANCE  OF  DIELECTRIC 
COATINGS  ON  SPHERICAL  RADIATORS. 


The  quantities  tabulated  are  the  ratios  of  the  expansion  coefficients  for 
the  field  of  a  coated  spherical  radiator  to  the  corresponding  coefficients 
for  the  field  of  an  uncoated  spherical  radiator  which  has  the  same 
excitation..  The  ratios  of  the  real  parts  and  of  the  imaginary  parts  are 
listed  separately.  The  parameters  are  as  follows:  sphere  radius  l.k 
wavelengths,  coating  thickness  .1  wavelength,  dielectric  constant  2.25. 
Data  and  output  tapes  are  on  file  at  the  U.  of  I.  Antenna  Laboratory  for 
other  parameters  as  follows: 

^   =      2.25,  t/A  -  .1,  -rA   :  .5(.l)l-6,  r/b   :  1.6(. 8)4.8; 

£T  -  +.o  ,  */*  =  ,1,  -T/h  :   .5(.l)i.6  ; 

£v  =  6.25,  V*  =  1.6,  */*  :  .l(.05)-3  . 


order 

FUcr 

IvnCw 

n 

R*o~ 

I^a^ 

1 

+15028  +01 

+11 5^1  +01 

2 

+11765  +01 

+178+3  +01 

3 

-1^931  +01 

+12138  +01 

k 

+1301+  +01 

+8ll40  +00 

5 

+10103  +01 

+17359  +01 

6 

+17206  +00 

+11925  +01 

7 

+27870  +01 

+86526  +00 

8 

+18860  +01 

++1378  +00 

9 

+30099  +01 

+2l8l4  -01 

10 

+1788+  +02 

+12306  +01 

11 

+1896^  +02 

+330+6  +01 

12 

+11221  +02 

+25209  +01 

13 

+91727  +01 

+219^6  +01 

14 

+85991  +01 

+2035^  +01 

15 

+85879  +01 

+19iK)3  +01 

8o 


TABLE    7  .     EIGENVALUES   OF  THE   SPHEROIDAL  WAVE 
Differential   equation   (1-v2)    dV2/dv2  +     ((piv)2  -  k    ) 


Pi  =5 


kjW 


odd 


1 

.21401687 

3 

.93590450 

5 

1 

.70632729 

7 

2 

.74589040 

.9 

4 

.10404858 

11 

5 

.78287310 

13 

7 

.78212780 

15 

10 

.10163380 

17 

12 

.74129183 

19 

15 

.70104605 

21 

18 

.98086382 

23 

22 

.58072513 

25 

26 

.50061724 

27 

30 

.74053164 

29 

35 

.30046262 

31 

40 

.18040621 

33 

35 

37 

39 

41 

43 

45 

47 

49 

51 

53 

Pi  =»  8 

HjAPi)2 

odd 

.12968845 

.58438654 

.98196932 

1.39873600 

1.92299702 

2.57428143 

3.35241653 

4.25662492 

5.28647244 

6.44172000 

7.72222957 

9.12791755 

10.65873113 

12.31463555 

14.09560733 

16.00163007 

18.03269202 

20.18878471 

22.46990174 

24.87603848 

27.40719101 

30.06335668 

32.84453316 

35.75071871 

38.78191194 

41.93811159 

45.21931685 


FUNCTION 

V  =   0      V(-l)    =   0 
Pi  =  12 

k/pi)2 
odd 

.08528034 

.39771381 

.67784891 

.92487780 

.1.16537067 

1.44669764 

1.78603032 

2.18356789 

2.63830632 

3.14963127 

3.71718528 

4.34075224 

5.02019593 

5.75542732 

6.54638607 

7.39303017 

8.29532979 

9.25326322 

10.26681421 

11.33597068 

12.46072310 

13.64106440 

14.87698883 

16.16849184 

17.51556974 

18.91821961 

20.37643922 


81 


(3i  =  5 


kjXpi)' 


odd 


55 

57 

59 

61 

63 

65 

67 

69 

71 

73 

75 

77 

79 

81 

83 

85 

87 

89 

91 

93 

95 

97 

99 

even 

2 

.58571819 

4 

1 

.29687773 

6 

2 

.18687186 

8 

3 

.38487933 

10 

4 

.90339175 

TABLE  7.  (CONTINUED) 

Pi  =8 

kj/(0i)2 

odd 

48.62552680 
52.15674084 
55.81295838 


even 

.36532339 

.78777844 

1.18065149 

1.64517729 

2.23271934 


Pi  =12 
kj/OiV 
odd 


21.89022638 
23.45957953 
25.08449723 
26.76497837 
28.50102197 
30.29262727 
32.13979344 
34.04251998 
36.00080629 
38.01465187 
40.08405647 
42.20901963 
44.38954113 
46.62562063 
48.91725799. 
51.26445283 
53.66720516 
56.12551471 
58.63938135 
61.20880496 
63.83378541 
66.51432263 
69.25041649 
even 

.24526082 

.54209080 

.80489795 

1.04264567 

1.29907596 
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Pi  =5 


j 

< 

k./([3i) 
sven 

12 

6 

.74246023 

14 

8 

.90185691 

16 

11 

.38144804 

18 

14 

.18115935 

20 

17 

.30094848 

22 

20 

.74078997 

24 

24 

.50066795 

26 

28 

.58057205 

28 

32 

.98049534 

30 

37 

.70043302 

32 

42 

74038172 

34 

36 

38 

40 

42 

44 

46 

48 

50 

52 

54 

56 

58 

60 

62 

64 

66 

68 

70 

TABLE  7.      (CONTINUED) 

PI  =  8 

kj/'((Bi)2 

even 

2.94755097 

3.78879486 

4.75586192 

5.84843150 

7.06632312 

8.40943005 

9.87768612 

11.47104867 

13.18948916. 

15.03298816 

17.00153171 

19.09511000 

21.31371552 

231.65734288 

26.12598793 

28.71964731 

31.43831867 

34.28199994 

37.25068942 

40 . 34438605 

43.56308857 

46.90679620 

50.37550829 

53.96922417 

57.68794341 


PI  =12 

^/(Pi)2 
even 

1.60901445 

1.97759588 

2.40383322 

2.88692212 

3.42639552 

4.02197700 

4.67349590 

5 . 38084249 

6.14394374 

6.96274954 

7.83722457 

8.76734344 

9.75308732 

10.79444239 

11.89139790 

13.04394555 

14.25207898 

15.51579320 

16.83508408 

18.20994839 

19.64038338 

21.12638702 

22.66795726 

24.26509285 

25.91779241 

27.62605496 

29 . 38987949 

31.20926533 

33.08421174 

35.01471814 
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Si   ■  5 


kjASi)' 


TABLE  7.      (CONTINUED) 
Si  =  8 


kj/(pi)' 


Si  =  12 


kj/CPi) 


72 

74 
76 
78 
80 
82 
84 
86 
88 
90 
92 
94 
96 
98 
100 


37.00078420 

39.04240935 

41.13959327 

43.29233563 

45.50063613 

47.76449461 

50.08391077 

52.45888438 

54.88941528 

57.37550344 

59.91714850 

62.51435054 

65.16710946 

67.87542497 

70.63929714 


FIGURE  13  RADIATION  PATTERNS,  Ha  VS.  0>OF  SPHERICAL  ANTENNAS  WITH  DIFFERENT 

CONDUCTING  SPHERE  RADII  XL         EXCITED  BY  EQUATORIAL  SLOT.  PATTERNS 
ON  LEFT  ARE  WT^h  no  COATING.   PATTERNS  ON  RIGHT  ARE  WITH  COATING, 


THICKNESS  .05  * 


£r  =  5 


f\ 


±L  z  1.61 


85 


=  2.41 


45°  90"  |35°  180°     ©  0 

FIGUEE  Ik    CONTINUATION  OF  FIGURE  13 
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RADIATION  PATTERNS, HW>  VS,  6  OF  SPHERICAL  ANTENNAS  WITH  DIFFERENT  CON- 
DUCTING SPHERE  RADII ,  rl3  EXCITED  BY  SLOT  AT  l60°  LATITUDE.   PATTERNS 
ON  LEFT  ARE  WITH  NO  COATING, 
THICKNESS  .1\,   e  =  2.25 


PATTERNS  ON  RIGHT  ARE  WITH  COATING, 


ii-  -.2.01 


87 


Sl  =  2.81 


0°    45° 

FIGURE  16 


90°     135°   180°  G       0*     45°    90     135°    180°   6 
RADIATION  PATTERNS,  Ha  VSe>  OF  SPHERICAL  ANTENNAS  WITH  DIFFERENT  CON- 
DUCTING SPHERE  RADIi;  T± 3  EXCITED  BY  SLOT  AT  l60°  LATITUDE.   PATTERNS 
ON  LEFT  ARE  WITH  NO  COATING,  PATTERNS  ON  RIGHT  ARE  WITH  COATING, 
THICKNESS  .IX,   €-2.25 


FIGURE  17  RADIATION  PATTERNS, H*  VS  e^OF  COATED  SPHERICAL  ANTENNA,-  =  1-51, 

COATING  THICKNESS  .1\,  FOR  DIFFERENT  VALUES  OF  DIELECTRIC  CONSTANT. 


SLOT  AT  160°. 


89 


RADIATION  PATTERNS, H^  VS  9,  OF  COATED  SPHERICAL  ANTENNA,  -  =  1.51> 
COATING  THICKNESS  .IX,  FOR  DIFFERENT  VALUES  OF  DIELECTRIC  CONSTANT. 
SLOT  AT  160°. 


90 


85*    WK 


FIGURE  19  RADIATION  PATTERNS,  H^  VS .  6,  OF  COATED  SPHERICAL  ANTENNAS,  v/X  =   1.51, 
DIELECTRIC  CONSTANT  6=  2.56,  FOR  DIFFERENT  COATING  THICKNESSES. 
SLOT  AT  160°. 


FIGURE  20  RADIATION  PATTERNS, %  VS.  e ^ OF  COATED  SPHERICAL  ANTENNAS. 
r/\  =  1.51,  DIELECTRIC  CONSTANT  €  =  2.56,  FOR  DIFFERENT 
COATING  THICKNESSES.   SLOT  AT  l60°. 
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e--o 


e--o 


Slot     '    a*    cen4< 


Uncoaled      a=  I.  00OO 


Coa+ed     a,--   l.ooooi 


FIGURE  21  PROLATE  SPHEROID  RADIATION  PATTERNS  (ONE  HALF  OF  SYMMETRIC 
PATTERN)  6  *  =   3,  £*.=  J5  i  H<>  «.  © 


nocoohng 

U^/07  7 


coded    u,H077 
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slot  of 

v0=-cos£o< 


PROLATE  SPHEROID  RADIATION  PATTERNS  p0  I 
DIFFERENT  SLOT  PATTERNS,  IUVS.  e 


no     coating 

U.  ,=    1.077 


coated 


u ,-  \   07  7 


Sk 


u,= 


1. 1 


FIGURE  23  PROLATE  SPHEROID  RADIATION  PATTERNS  p  i  =  5,  € 
DIFFERENT  SLOT  PATTERNS,  H.  VS.  6 


&,   FOR 


no  coating 
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!     o 


FIGURE  2k     PROIATE  SPHEROIDAL  RADIATION  PATTERNS,    PQ    I 


H,  VS    0 


no  coating 
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Vo-'COSdO0 


135°     180°  e 


FIGURE  25     PBOIATE  SPHEROIDAL  RADIATION  PATTERNS,    pQ  I    =  Q,   Hx  VS. 6 


no    coating 

TA^s- 1.077 


Mote: 

concJuctoY  si*e.  tkat  of  Fiq.  aa 


V.-  "  COS  30 


*5*    90°     (35°    160'     ^ 

FIGURE  26  PROLATE  SPHEROIDAL  RADIATION  PATTERNS  3A=5"j  12,  H^  VS.  e 


no  coating 


\      Vo  =  -C0S45° 


j/ot 
at 
Vo^-coseo° 
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coa"ted 
Note: 


condacxot  siie  thai  ©?-  Fit.  jl3L 
11,.=  1.077 


^*   S06    A35-*   /5b 


FIGURE  27  PROLATE  SPHEROIDAL  RADIATION  PATTERNS  g£-  5512,  Ha  VS.© 
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100 
NOTE:. 

In  Figures  29  to  63,  negative  numbers  are  represented  on  the  semi- 
log scales  "by  the  following  conventions.  (The  lines  joining  points  are 
merely  aids  for  the  eye) 

Real  Part 

0    positive  value 
X    negative  value 

joins  positive  values 

.  .  ,  joins  negative  values  or  positive  and  negative  values 

Imaginary  Part 

•     positive  value 
J7]    negative  value 

joins  positive  values 

I  I I  J  I  \   I  J  I  .  joins  negative  values  or  positive  and  negative  values 
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—  REAL  PART 

—  IMAGINARY  PART 
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FIGURE  29.[H/(2jtr/\)]"     AS  A  FUNCTION  OF  n 
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REAL  PART 
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FIGURE  30.[H/(2rtr/\)]"1  AS  A  FUNCTION  OF  n 
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REAL  PART 
■IMAGINARY  PART 
(r/\  =  1.0) 


13  5  7  3  1/ 

FIGURE  31-[H/(2nr/\)]"1  AS  A  FUNCTION  OF  n  . 
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REAL  PART 

IMAGINARY  PART 

(r/\  =  1.2) 


FIGURE  32.[H'(2jtr/\)]   AS  A  FUNCTION  OF  n 
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REAL  PART 
IMAGINARY  PART 
(r/\  -  1A) 


FIGURE  33. [E'{2xr/\)]        AS  A  FUNCTION  OF  n 
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REAL  PART 
IMAGINARY  PART 
(r/X  =  1,8) 


n 


5      A  9        13        i"? 
FIGURE  35,  [H,(2jtr/\)]_  AS  A  FUNCTION  OF  n. 
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REAL  PART 
IMAGINARY  PART 
(r/\  =  2.0) 


FIGURE  36.  [H/(2itr/\)]   AS  A  FUNCTION  OF  n 


109 


1 


REAL  PART 
IMAGINARY  PART 
(r/\  =  2.2) 


1 

o 

1 

\     I 

>\    / 

( 

i  JF  ' 

1     • 

A  i 

\  / 

1  J 

\     l\ 

/  • 

V 



',  ft  1 

\    1 

•  i  / 

V 

ft 

/ 

^ 

/(' 

» 1      / 

V     1/ 

/ 

I 

/ 

1 

J  1 ' 

T  T 1 

J  y 

/   , 

1 
1    1 

/ 
/ 

I 

1  I' 

\  /I 

/■ 

1 

1 

1 

\ 

Ji      li 

r 

1 1 

*  1 

/ 

1      \ 

\ 

i 

\ 
\ 

1 

■!■ 

( 

/ 

(     ; 

: 

i 

/ 

\ 

< 

\ 

,1 

i  1 

it 

u 

3 

' 

i 
I 
l 
i 

I 

i 

1 

i 

1 
1 

' 

1 
1 

_ 

1 

1 

\ 

> 

\ 
\ 
\ 

01 

\ 
\ 
\ 
\ 
\ 

\ — 

r\ 


3    5    7    9    M    13    15    \7   19 
FIGURE  37.  [H/(2jtr/\)]_1  AS  A  FUNCTION  OF  n  . 
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REAL  PART 
IMAGINARY  PART 
(r/\  =  2.1+) 


5  7         9         II         13        15        II         19      21 

FIGURE  38.  [H/(2itr/\)]"1  AS  A  FUNCTION  OF  n  . 
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REAL  PART 
■IMAGINARY  PART 
(r/\  =  3.2) 


3    5 


7    9    II    12    15    n    f9   21    It)   2.5  2( 
FIGURE  39.  [H/(2nr/\)]"1  AS  A  FUNCTION  OF  n. 
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6        7       8        9       10      II      n- 

FIGURE  kO     RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE  COATED 
SPHERICAL,  RADIATOR  TO  THE  COEFFICIENTS  a  IN  THE  EXPANSION  FOR 
THE  UNCOATED  SPHERICAL  RADIATOR  WITH  THE  SAME  EXCITATION 
(g  =  r/\)   COATING  THICKNESS  .1\.  €  =  2,25 
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FIGURE  41  RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE  COATED 
SFHEROECAL  RADIATOR  TO  THE  COEFFICIENTS  a   IN  THE  EXPANSION  FOR 
THE  UNCCATED  SPHERICAL  RADIATOR  WITH  THE  SAME  EXCITATION. 
(g  =  r/\).      COATING  THICKNESS  .IX-   €  =  2.25 
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FIGURE  k-2     RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE  COATED 
SPHERICAL  RADIATOR  TO  THE  COEFFICIENTS  a   IN  THE  EXPANSION  FOR 
THE  UNCOATED  SPHERICAL  RADIATOR  WITH  THE  SAME  EXCITATION. 
(g  =  r/x).   COATING  THICKNESS  .IX.   6=  2.25 


FIGURE  kj>     RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE  COATED 
SPHERdCAL  RADIATOR  TO  THE  COEFFICIENTS  a   IN  THE  EXPANSION  FOR 
THE  UNCOATED  SPHERICAL  RADIATOR  WITH  THE  SAME  EXCITATION. 
(g  =  r/\).   COATING  THICKNESS  .1\.   €=  2.25 
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RATIO  OF  COEFFICIENTS,    c    ,    IN  THE  FIELD  EXPANSION  FOR  THE  COATED 

SPHERICAL  RADIATOR  TO  THE  COEFFICIENTS  a      IN  THE  EXPANSION  FOR 
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THE  UNCOATED  SPHERICAL  RADIATOR  WITH  THE  SAME  EXCITATION. 
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FIGURE  58  RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE  COATED  SPHE- 
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FIGURE  59  RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE 
COATED  SPHERICAL  RADIATOR  TO  THE  COEFFICIENTS,  a  ,  IN 
THE  EXPANSION  FOR  THE  UNCOATED  SPHERICAL  RADIATOR1  WITH 
THE  SAME  EXCITATION.   €=  2.25,  £  =  1*6.   THICKNESS  t/\ 
VARIED  IN  SUCCESSIVE  GRAPHS 
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FIGURE  60  RATIO  OF  COEFFICIENTS,  c  ,  IN  THE  FIELD  EXPANSION  FOR  THE  COATED 
SPHERICAL  RADIATOR  TO  THE  COEFFICIENTS,  a  ,  IN  THE  EXPANSION 
FOR  THE  UNCOATED  SPHERICAL  RADIATOR  WITH  THE  SAME  EXCITATION. 
€  =  2.25,  £  _  1.6.   THICKNESS  t/\  VARIED  IN  SUCCESSIVE  GRAPHS 
\  ~ 
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APPENDIX  A 

Determination  of  the  Eigenvalues  of  Unfamiliar  Eigenfunctions 

The  purpose  of  this  appendix  is  to  outline  a  general  method  for 
the  calculation  of  the  eigenvalues  which  are  associated  with  a  second 
order  differential  equation  with  boundary  conditions,  in  the  event 
that  the  functions  which  satisfy  the  differential  equations  are  unknown. 
This  problem  is  different  from  the  problem  of  eigenvalue  determination 
as  it  is  usually  encountered  in  engineering  problems  in  that  the 
functions  which  represent  the  solutions  are  known  only  in  the  form  of 
some  kind  of  infinite  series,  and  the  dilemma  is  that  the  coefficients 
in  the  infinite  series  depend  on  the  eigenvalues  (which  are  as  yet  unknown) 

The  method  discussed  here  is  essentially  a  generalization  of  the 
nfethod  described  in  Section  2.2.3  for  the  calculation  of  the  eigenvalues 
of  the  spheroidal  functions,   (it  is  assumed  that  the  reader  is  familiar 
with  the  idea  that  orthogonal  functions  are  generated  by  differential 
equations.  This  is  discussed  briefly  by  Churchill  and  thoroughly  by 
Ince  . )  The  method  is  probably  well  known  to  mathematicians,  but  it  is 
not  presented  specifically  in  any  of  the  books  into  which  the  engineer 
ordinarily  looks  for  help  in  solving  his  problems . 

A.l     Consider  a  general  second  order  linear  differential  equation  of 
the  Sturm-Liouville  type 


p(x)y  +  q(x)y  +  (r(x)  +  kks(x))y  = 


(Al) 
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with  the  boundary  conditions 

^y(a)  +  p.,/ (a)  =  0 

(A2) 
a2y(h)   +  fi2y'(-b)   =  0 

The  problem  is  to  find  the  eigenvalues,  k,    ,  in  the  event  that  the 

functions  y   are  unfamiliar.  To  find  these,  expand  the  functions  y 
k  k 

in  a  series  of  known  functions  which  form  a  complete  orthogonal  set  and 
satisfy  the  equation 

P^x)^"  +  q1(x)^n/  +  (^(x)  +  \nSl(x))</>n=0  (A5) 

and  identical  boundary  conditions  (A2) .   (That  is,  the  eigenvalues  of 

(f)      are  known  or  can  be  readily  found.)  Thus,  assuming  that  term  by 
n 

term  differentiation  is  possible,  the  series 


^  ■  I  Ckn* 
n=l 

is  substituted  into  (Al)  with  the  result 


L    Vn  W) 

n=l 


*<x>  Y   Ckn<  +  *«  I  CknV  +  <r<x)  "  V«>  }_  ClA  "  °  .  (A5) 
n  n  n 

To  generate  a  system  of  linear  equations  having  the  form 


y 


C.   L   =0 
kn  nm 


multiply  through  by  d>     and  integrate  over  the  range.  The  result, 
m 

after  interchanging  the  order  of  summation  and  integration  and  grouping 

the  coefficients  of  c       ~>   is 
kn 


ltei 


i 


'kn 


P(x)^n<^mdx  +    /       q^ty^dx  +        rUW^dx+kJ     s^^dx 
Ja  Ja  ^  (A6) 


=0 


The  quantity  in  brackets  is  an  array  of  numbers  which  will  be  designated 

L       .  The  eigenvalues  are  those  numbers  which  make  the  determinant 
of  this  array  equal  to  zero,  since  this  is  the  condition  for  consistency 
of  the  set  of  equations.  As  is  shown  in  Section  A. 2,  the  values  obtained 
from  a  finite  size  determinant  will  be  good  approximations  since  the 
representation  has  stationary  properties. 

The  determination  of  the  elements  of  L      can  be  simplified,  and 

nm  *  ' 

the  rate  of  convergence  improved,  by  the  judicious  choice  of  the  functions 
0      .To  show  this,  multiply  each  member  of  the  set  of  Eqs.  A3  by 


the  corresponding  member  of  the  set  of  constants,  C. 
equations  together  to  form 


kn 


,  and  add  these 


n  n  n 

Now  subtract  this  equation  from  (A5),  multiply  successively  by  the  set 


(A7) 


of  functions,  <$>  ,  and  integrate  over  the  range  of  orthogonality  to 

obtain  the  result 
I—  „b 


dx 


Y_   Ckn     |   <*<*>  "  PiM^V*  +  /   ((l(x)  "  ^Vn^W 
n        J  a  <J   a 


(r(x)  -  ri(x))0n0mte  -  yi^  +  K  iW^V 


in  which  N      is  the  normalization  integral  of  the  functions,    0 


(A8) 


that  is 


Un 


N  4 


nm   w  nm 


s(X)0n^dx 


The  elements  of  L     in  the  brackets  of  (A8)  may  be  easier  to  evaluate 
than  those  of  (A6) .   In  particular,  Eq.  A8  displays  formally  the 
(intuitively  obvious)  conditions  which  make  (f>         a  favorable  set  of 
functions,  namely,  they  should  satisfy  a  differential  equation  which 
is  as  nearly  like  the  given  equation  as  possible.  The  best  possible 
condition  is  that  all  except  one  of  the  following  hold: 

p(x)  =  px(x)  ,  q(x)  =  qx(x)  ,  r(x)  =  r^x)  ,  s(x)  =  s^x)  . 


Otherwise,   0    should  be  selected  so  as  to  make  the  integrations  as 
simple  as  possible. 

We  note  also  from  (A8)  the  condition  for  symmetry  of  the  matrix 

|l   =  L    ): 

nm    mn 


J 


(p-p^/c^dx  +      (cL-q1)^i/0ndx. 
J  a 

(A9) 


Thus,    if     p-p  and        q  =  q         >   there  is  no  question  of  the  symmetry 

of  the  matrix,      t  .     An  integration  by  parts  permits    (A9)   to  be 

nm 

replaced  by  the  condition 


r 


L' 


£_(p-P1)-(q-q1) 

dx 


<Mr>  =  [^(P-Pi)^^! 


4>  4  ' 

TxFm 


dx, 


(A10) 


This  is  a  sufficient,  but  perhaps  not  a  necessary  condition. 


Ik2 

There  are  a  few  more  comments  which  can  he  made  ahout  (A6)  and 

(A8).   If  the  values  of  the  elements  so  calculated  were  placed  in  an 

infinite  determinant,  and  the  condition  I  L   1=0     enforced,  then 

J   nm 

exact  eigenvalues  would  he  obtained .  But  the  vital  question  is  whether 

finite  determinants  can  he  employed  to  obtain  usefully  accurate  values. 

The  use  of  such  a  finite  determinant  can  he  regarded  as  equivalent  to 

the  replacement  of  the  infinite  series  (A^-)  hy  a  finite  sum  of  N  terms 

as  follows : 

N 


yk  =  ^   Ckn*n   ' 


n=l 

The  practical  questions  are:  Can  the  eigenvalues,    k_       ,  he 

obtained  to  arbitrary  accuracy  with  a  finite  N  ,  and  if  so,  how  large 

must  n    be  to  obtain  a  specified  accuracy?  The  answer,  of  course, 

is  inseparably  connected  to  the  matter  of  how  closely  the  functions 

y,     and   0,    resemble  each  other.  In  any  case,  however,  as  is 

shown  in  Section  A .2,  (A6)  and  (A8)  are  results  which  have  stationary 

properties  (that  is,  they  are  the  same  results  as  obtained  by  the 

application  of  a  variational  procedure).  Thus,  the  end  result  is 

N 
insensitive  to  small  differences  in  the  quantities  y,     and  2__  Ckn^n 

n=l 
Moreover,  the  assumption  of  uniform  convergence  is  necessary  in  order  to 

justify  the  term  by  term  differentiation  to  arrive  at  (A5)»  Such  uniform 

convergence  means  exactly  that 
N 


I'k"  ZC*n'>'n!  -' 


€ 

n=l 


1*3 


independent  of   x   ,  for  some  finite   N   .  This  implies  that 

C   =  0     >       M>?  N       ,  so  that  the  accuracy  of  the    th 
km 

eigenvalue  can  be  determined  to  an  accuracy  set  by  the  value  of  € 
selected.  In  practice,  the  error  analysis  is  difficult,  so  that  with 
a  high  speed  computing  machine  available,  the  question  is  perhaps  best 
answered  by  repeated  trials  in  which  the  size  of  the  matrix  is  varied. 
It  is  interesting  that  in  the  case  of  the  spheroidal  functions  (Section 
2.2.3);  even  the  value   N  =  1   gives  a  correction  to  the  eigenvalue 
of  the  trial  function  P       (lowest  eigenvalue). 

The  idea  of  replacing  the  infinite  sum  (A^)  by  a  finite  sum  suggests 
one  other  technique  which  is  useful  in  practice  in  the  calculation  of  the 
higher  order  eigenvalues.  For  with  the  higher  order  functions,  a  good 
approximation  may  be  as  follows : 
N 


yk 


/  kn  n  , 


that  is,  the  lower  order  coefficients  may  also  be  approximately  zero. 
This  means  that  the  accuracy  of  the  higher  order  eigenvalues  can  be 
improved,  with  a  given  size  matrix,  by  taking  a  section  "out  of  the  middle" 
of  the  infinite  matrix  (that  is,  dropping  out  some  of  the  lower  order 
coefficients  as  well  as  some  of  the  higher  order) .   (Of  course,  the 
coefficients  for  all  orders  of  the  functions  satisfy  the  same  set  of 
(infinite)  equations,  so  the  solution  of  a  single  matrix  gives  an 
approximation  for  as  many  eigenvalues  as  coefficients  included  — 
however,  the  accuracy  of  the  individual  eigenvalues  so  determined 


Ikk 

depends  primarily  on  how  well  the  finite  sum  for  each  order  approximates 
the  function  of  that  order. 

A. 2     It  will  now  he  shown  that  a  variational  method  leads  to  Eq.  A6 
and  hence  to  the  same  determinant  as  the  one  ohtained  earlier.  Consider 
the  differential  equation  in  operational  form 

Ly  +  *y  =  0    • 

Multiply  this  equation  by  a  trial  function,   *\>    ,   and  integrate  over 
the  complete  range: 

/  if  Lydx  +  \  j  \\>   ydx  =  0   • 

Now  perform  a  first  variation 

JS^  Lydx  +  S   \ft  ydx  +  \  W  ydx 

and  set  £  \   =  0    to  obtain  the  result 
J  (Ly  +  \y)j>  dx  =  0   . 

In  the  present  problem,  Ly  +  \y   is  of  course  just  the  right  hand 
side  of  Eq.  Al,  which  implies 

J  p(x)y"  +  q(x)y'  +  (r(x)  +  kks(x)y)^  dx  =  0 

Now  assume    y  =  )       C  <$>         .   take^   =  )     C  <t>  (so/     =  )  /  C  0  ) 

*      /_ n  n  / mm  /  °    m  m 

n  m  m 

and  find 

p(x)  Y_*X    +  4(X)   L  CA'    +   <r«  +  VW][CaJX[iC.e 


1*5 


orJjCm[       (p(x)  ^'CA*  +  q(x)  Z  *'  +   (r(x)   +  kks(x))ZCnt)(^rH=0  ' 
m          O a                  n  n  n 

But  the  variations  in  C  >  £   C       ,  must  be  arbitrary,    so 
b 

Cp(x)^  Cn0/  +  q(x)2_  Cn0/    +   (r(x)   +  kks(x))^cn0n>mdx  =  0 


L n 


p(x)<k"0m  +  <l(x)0n>m  +  (r(x)  +  VCxJ^^  dx  =  0  , 


n   *^a 
which  is  the  result  stated  in  Eq..  A6.  Thus  we  have  demonstrated  that 
(A6)  has  stationary  properties. 


ll+6 
APPENDIX  B 

Orthogonality  Properties  of  the  Vector  Functions  of  Chapter  3 

In  the  problem  of  the  spherical  antenna  excited  by  an  arbitrary 
slot  configuration ,    it  is  convenient  to  represent  the  transverse  electric 
field,   E.-r   ,  by  means  of  an  expansion  in  orthogonal  functions.   In  order 
that  this  can  be  done^  the  following  set  of  othogonality  relations  is 
required: 


/e\      .  e\  dS  =  L  .  1  iV 
le    je      ij  J  it 


H/    2 

p  dS 
le 


I  j 

r=rn 


e.      .   e.  dS  =   d  , .    J  e. 

J  ie         je  ij  J     i< 


eio   •  V3  =  <*  i  J   el< 

r=rl 

I 


(Bl) 


f  C-^^ijC*3  <b2> 

j        e"ie  •  SJo«  "  °  <B5> 


e;    .  ev  dS  =  0  (B^) 


(B5) 


2ds  (B6) 


eie   •  •to*  =  °  (B7) 


1^+7 

in  which 

e 'e  =  V  T  « 

1 0  1 '-"  , 

e.e  =  v  T.e  X  r   . 

The  relations  (Bj5)  and  (B7)>  which  involve  products  of  the  even  and 
odd  vector  functions,  follow  almost  immediately  since  the  integrals 
will  all  involve  the  product  of  a  sine  function  and  a  cpsine  function 
integrated  over  the  range  0    to  2n 

The  relationship  (B^)  can  be  shown  in  general  by  employing  the 

divergence  theorem  and  some  well  known  vector  identities  (listed  for 

example  in  Stratton,  p.60^)  as  follows: 
J  V*  Adv  =  J  A*nds 

put       A  =  T.(r  xv  T.)  .  Then 

■*■      J 

V  •  A  =VT.  •  (r  xVT.)  and 


J  A  .  ndS  =J  T.ffixVT.)  • 


rdS 


VT.  s  (rxyT  )dS  =   e.  .  e.  dS  =  0 

The  relation  (Bl)  can  be  demonstrated  as  follows^-  The  terms  ofyT. 
m  a  ie 

are  Q   d  n   cos  m<f)  -<j)m         P  m  sin  m  cp      5         so  that  the  product 

r  d©  r  sin<9  n 

of  V  T .   »yT,         is  composed  of  such  terms  involving  all  possible 

combinations  of  m   and   n   .   It  is  clear  therefore  that  unless  the 

index   m  is  the  same  in  bothvT,      and-  VT-    >   "the  <r  -integration 

will  give  zero.  Thus  the  product  will  be  of  the  form 

1  dP  m   dPnm      2  ,    2     ^  m  „  m  .  2   1 
~o-   n     %         cos  mp  +  m P   P.   sin  m  (f> 

r    i&  &&  2  ".    2   n   X 

r  sin  © 


XkQ 


After  the   $>      -integration,  the  results  are  integrals  of  the  type 


,  dP  ™  dP„       2  _  m  „  mx   .  ^  a 

(   n     I         +     m  P   P0  )  sin#d0 

15"  TT    "5  n   * 

sin6> 


This  integral  is  zero  unless   n=  I  >  and  constant  otherwise 

(Stratton,  p.  ^17). 

The  demonstration  of  (B2)  is  completely  analogous.  Furthermore, 

the  results  (B5)  and  (B6)  follow  from  the  same  type  of  argument  since 

the  operation  £  x  ^p  e        (which  is  the  operation  required  to 
i  o 

transform  £  e       into  i  e       )  merely  interchanges  the  O 
i°  j© 

and  0      -components. 


ERRATA 


Page 
iv        In  title  of  Fig.  22,  23,  replace   €  =  JZ   by  € 

iv        In  title  of  Fig.  28,  replace  —~. ^ —    by 


kl 


48 


6o 

69 

93, 

9k 

112- 

130 

f '(p  f  Ih'O  f 

n  r'  I  n  r  I 

iv-vlii     Add  one  to  the  page  number  for  j>ig.  29  and- beyond.  Page  100  is 
note  explaining  the  symbols  in  Fig.  29  and  those  beyond. 

16        Line  8  insert  (Columbia  University  Press  19^+7)  after  Tables  of 
Spherical  Bessel  Functions. 

18        Next  to  last  line;  replace  (33)  by  (23). 

A 

23        All  quantities  Z  should  read  Z  . 

2k  Last  line;  replace  "the  lower  order  ratio"  by  "the  next  lower 

order  ratio". 

a' 

26        Subscript  n  left  off  N  in  two  equations. 

kl  Line  21;  replace  appraicably  by  appreciably. 

Line  22 j  replace  magnitude  by  magnitudes. 


Line  17;  replace   /  d,    2n(n+l)  by 


lkn  2n<n+1) 


2n+l 

10  10 

Center  block;  replace    )   i   P  by    )   cL   P 

n=l  n=l 

Third  box  from  top,  on  right,  replace     V7  d   2n(n+l)    by 
V-2  Ls  2n+l 

n 
Line  3;  replace  (62)  by  (63). 

/ 
Equation  89  replace  d   g  by  d  %      . 

First  line  after  equation  89.   Insert  e.  at  the  end  of  the  line. 

/  a 
Fig.  22,  23.  Replace    €^     =    ■=     by  €^ 

Ordinate  of  Figures  kO-^d   should  all  be  lower  case,  /  n 


ERRATA  (continued) 

]>3       Line  9.   Replace  P  "by  P- 

1^J+       Next  to  last  line;  replace  6   =    >   6  C  0   by 6^ =  /  6C  0 

m  m 

1^5       First  equation;  replace  q(x)  }       <p/    by  q(x)y   C  0/ 
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